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1.0  Introduction 


The  physical  processes  that  occur  in  the  boundary  layer  along  a  streamlined  body 
at  incidence  are  very  complex  and  far  from  being  completely  understood.  These  include 
such  phenomena  as  transition  from  laminar  to  turbulent  flow,  the  presence  of  a  three- 
dimensional  laminar  or  turbulent  boundary  layer  with  strong  crossflow,  three-dimensional 
laminar  or  turbulent  boundary  layer  separation,  and  the  formation  of  longitudinal  vortices 
extending  into  the  wake. 

The  successful  understanding  and  modeling  of  these  physical  processes  can  pro'-? 
fruitful  in  many  applications.  The  ability  to  control  transition  and  separation  can  aid  in 
decreasing  the  aerodynamic  drag  on  a  body.  Improvements  in  the  maneuverability  of 
aircraft,  missiles,  and  submarines  can  result  from  the  understanding  of  the  processes 
involved  in  three-dimensional  turbulent  boundary  layers  and  the  formation  of  longitudinal 
vortices. 

Solutions  of  the  Navier-Stokes  equations  car.  yield  accurate  predictions  for 
complex  flows,  including  separation;  however,  such  solutions  would  require  tremendous 
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computer  resources  and  their  accuracy  is  dependent  upon  results  from  turbulence  models. 
In  the  past,  attempts  to  develop  three  dimensional  boundary  layer  models  including 
separation  and  trans:tion  have  not  met  with  great  success  and  the  extension  of  two- 
dimensionM  models  to  three-dimensional  situations  has  not  provided  consistently  useful 
results.  For  the  time  being,  experimental  data  must  be  used  as  a  basis  for  understanding 
these  complex  flows. 

The  great  interest  in  transition  and  separation  phenomena,  and  in  flow 
characteristics  in  the  presence  of  vortical  flow,  is  shown  by  the  amount  of  experimental 
research  that  has  been  dedicated  to  investigating  the  processes  that  occur  in  these  flows. 
A  number  of  experimentalists  have  investigated  flow  around  generic  axisymmetric  bodies, 
such  as  a  6:1  prolate  spheroid,  because  the  flow  field  around  such  a  body  is  made  up  of 
the  same  features  that  are  of  interest  when  considering  the  flow  field  around  an  aircraft 
fuselage,  a  missile,  or  a  submarine. 

Research  conducted  at  DFVLR  in  West  Germany  on  a  2.4  meter  long  6:1  prolate 
spheroid  has  obtained  a  great  amount  of  experimental  data.  Wall  shear  stress 
measurements,  surface  oil  flow  visualization,  surface  pressure  measurements,  and  mean 
velocity  field  measurements  have  yielded  information  about  boundary'  layer  transition  and 
separation  (Kreplin,  et  al.,  1980.  1982;  Meier  and  Krephn,  1980;  Meier,  et  al.,  1981a, 
1981b.  1982,  1983a,  1983b).  Mean  boundary  layer  profile  measurements  are  presented 
in  Meier,  et  al.  (1984),  and  outer  mean  velocity  field  measurements  are  given  in 
Vollmers,  et  al.  (1985).  Results  of  some  preliminary  turbulence  measurements  obtained 
at  DFVLR  are  presented  in  Kreplin  and  Vollmers  (1981).  Table  1  shows  the  flow 
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conditions  for  the  various  measurements  conducted  at  DFVLR.  The  results  of  all  these 
experiments  are  used  as  a  test  case  by  other  researchers. 

Cooperative  research  conducted  between  DFVLR  and  ONERA  has  resulted  in 
detailed  investigations  of  turbulent  separation  on  an  oblate  spheroid,  including  surface 
flow  visualization,  surface  pressure  measurements  and  velocity  field  measurements 
(Chanetz.  1986:  Barberis  and  Chanetz.  1986:  Chanetz  and  Delery.  1988). 

A  Sino-Italian  research  program  investigating  the  three-dimensional  flow  around 
a  6:1  prolate  spheroid  has  yielded  results  for  flow'  visualization,  pressure  and  force 
measurements,  boundary  layer  measurements,  and  velocity  field  measurements  (Iuso.  ei 
a!..  1988). 

Besides  investigating  transition  and  separation  phenomena  around  an  axisymmetric 
body,  a  number  of  researchers  have  studied  the  interactions  of  a  vortex  and  a  turbulent 
boundary  layer.  Westphal.  el  al.  (1985)  investigated  such  interactions  in  the  presence  of 
a  streamwise  pressure  gradient.  Pauley  and  Eaton  (1988.  1989)  studied  the  effects  of 
vortex  pairs  embedded  in  a  turbulent  boundary  layer.  Other  work  has  been  don  by 
Bradshaw  and  his  co-workers  (Cutler  and  Bradshaw.  1986:  Shabaka.  Mehta  and 
Bradshaw  ,  1985;  Mehta  and  Bradshaw.  1988).  All  of  these  investigations  include  mean 
velocity  and  Reynolds  shear  stress  measurements  to  study  the  behavior  of  the  flows. 

This  research  is  part  of  an  overall  investigation  of  the  flow  field  around  a  6:1 
prolate  spheroid,  in  both  steady  and  unsteady  flows.  The  present  study  involves  the 
investigation  of  the  region  of  separation  on  the  leeward  surface  of  the  6:1  prolate 
spheroid  at  moderate  angles  of  attack  in  a  steady  flow.  Mean  velocity  field  measurements 
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in  the  near  wake  region  with  the  model  at  angles  of  attack  a  =10°  and  15°  and 
Re=1.3xl06  and  4.0xl06  allow  the  analysis  and  description  of  the  structure  of  the 
vortical  flow  field.  Mean  velocity  and  turbulence  measurements  in  the  boundary  layer 
at  a  =10°  and  Re  =  4.0 x  106  allow  the  analysis  and  description  of  the  three-dimensional 
turbulent  boundary  layer,  including  the  growth  of  the  boundary  layer  along  the  surface 
of  the  model,  and  any  crossflow  and  separation  phenomena  that  are  evident.  The 
turbulence  measurements  are  used  to  attempt  to  appraise  the  extension  of  two-dimensional 
turbulence  models  to  three-dimensional  flows.  The  choices  of  angle  of  attack  a  and  Re 
are  based  upon  the  results  of  the  mean  flow  field  measurements  and  the  surface  oil  flow 
measurements  which  indicate  that  at  these  conditions,  flow  phenomena  leading  to 
separation  occur  over  an  extended  region  along  the  surface  of  the  model. 

The  results  of  this  study  can  possibly  be  used  as  a  basis  for  testing  the  validity  of 
three-dimensional  turbulence  models  by  providing  experimental  data  for  comparison. 
More  importantly,  the  results  can  be  useful  for  comparison  when  future  measurements 
are  made  for  unsteady  flow  around  the  same  6:1  prolate  spheroid  model,  as  part  of  a 
continuing  investigation  of  transition  and  separation  processes. 
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2.0  Experimental  Apparatus  and  Instrumentation 

2.1  Wind  Tunnel 

All  the  measurements  were  conducted  in  the  Virginia  Polytechnic  Institute  and 
State  University  Stability  Wind  Tunnel.  This  wind  tunnel  is  a  continuous,  closed  test 
section  subsonic  w'ind  tunnel  with  interchangeable  round  and  square  cross-section  test 
sections.  Originally  constructed  in  1940  at  NASA  Langley  in  Hampton,  Virginia,  the 
tunnel  was  moved  to  Virginia  Tech  in  1958. 

As  shown  in  Figure  1,  this  wind  tunnel  is  a  closed-loop  tunnel  powered  by  a  600 
HP  D.C.  motor  driving  a  propeller  4.27  meters  in  diameter.  This  provides  a  maximum 
speed  of  67  mps  in  a  normal  six  feet  by  six  feet  test  section  configuration.  An  air- 
exchange  tower  located  downstream  of  the  fan  and  motor  assembly  provides  temperature 
stabilization.  Each  comer  of  the  wind  tunnel  circuit  has  streamlined  turning  vanes  spaced 
along  the  diagonal  of  the  comer  to  provide  uniform  flow.  Upstream  of  the  test  section, 
in  the  settling  chamber,  are  seven  turbulence  screens.  Betw-een  the  test  section  and  the 
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settling  chamber  is  a  nozzle  with  a  9:1  contraction  ratio.  Downstream  of  the  test  section 
is  a  diffuser  with  an  angle  of  3°.  Vortex  generators  located  at  the  beginning  of  the 
diffuser  aid  in  efficient  diffuser  operation. 

The  test  section  used  in  these  measurements  is  23  feet  long,  with  a  six  feet  by  six 
feet  square  cross  section.  An  air-tight  test  and  control  room  encloses  the  test  section  to 
minimize  air  leakage  into  the  test  section.  The  6:1  prolate  spheroid  model  was  mounted 
on  a  dummy  balance  and  sting  and  supported  in  the  test  section  by  a  six-feet  tall  strut. 
The  strut  and  sting  allowed  the  model  to  be  set  at  angles  of  attack  between  approximately 
-15°  and  +22°. 

2.2  6:1  Prolate  Spheroid  Model 

The  6:1  prolate  spheroid  model  used  in  these  experiments  is  1.37  meters  (54 
inches)  in  length  and  has  a  maximum  diameter  of  0.23  meters  (9  inches).  The  model, 
showm  in  Figure  2  and  3,  consists  of  a  forward  and  aft  section  of  identical  design,  with 
an  aluminum  internal  framework  of  four  bulkheads  and  eight  longerons  surrounded  by 
a  fiberglass  skin.  The  two  sections  are  joined  together  by  a  center  section,  an  aluminum 
cylinder  covered  by  fiberglass  skin,  that  can  accommodate  a  balance  and  sting. 
Interchangeable  aluminum  nose  and  tail  pieces  are  mounted  at  the  respective  ends  of  the 
forward  and  aft  sections.  The  tail  piece  allows  the  sting  and  balance  to  pass  through  the 
hollow’  interior  of  the  model  for  mating  with  the  center  section. 
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2.3  Tunnel  Traverse  and  Coordinate  System 

The  probes  used  in  all  the  measurements  were  mounted  on  a  traverse  unit  that 
allowed  movement  in  both  horizontal  and  vertical  directions.  These  directions  correspond 
to  Zj-  and  YT.  respectively,  as  defined  by  the  right-hand  orthogonal  tunnel  coordinate 
system,  shown  in  Figure  4.  used  throughout  the  experiments.  The  circumferential  angle 
<p  is  measured  with  respect  to  the  -Zp  axis  and  is  positive  in  the  direction  of  the  YT  axis. 
The  traverse  mechanism,  controlled  by  two  separate  manually-operated  control  units, 
allowed  a  minimum  movement  of  0.254  mm  (0.01  inches)  in  the  horizontal  (Zp)  direction 
and  0.1  mm  (0.0039  inches)  in  the  vertical  (YT)  direction.  Before  each  run,  the  probes 
were  placed  at  the  desired  location  in  the  XT  direction  by  changing  the  position  of 
adjustable  brackets  used  to  mount  the  probes  on  the  traverse. 


2.4  Three-Dimensional  Yawhead  Probe  Instrumentation 

•  A  United  Sensor  DC  125  three-dimensional  yawhead  probe  (United  Sensor 

Catalog)  w’as  used  to  obtain  mean  velocity  field  measurements  in  the  near  wake  region. 

This  probe  has  a  blunt  conical  head  with  five  pressure  ports,  one  port  in  the  center 

ft 

surrounded  by  four  ports  spaced  90°  apart.  The  magnitude  of  the  velocity  and  its  pitch 
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and  yaw  angles  are  found  by  measuring  the  pressure  at  all  five  ports.  The  probe 
calibration  used  (Sung,  1985)  allowed  the  determination  of  velocity  vectors  within  a  60° 
cone  from  the  axis  of  the  probe. 

The  probe  was  mounted  on  the  Y T-7T  traverse  described  in  Section  3.3.1,  and 
oriented  parallel  to  the  Xx  direction.  Any  probe  misalignment  with  the  XT  axis  was 
measured  using  a  cathetometer  and  a  digital  protractor,  and  corrections  were  made  during 
the  data  reduction.  The  position  of  the  probe  with  respect  to  the  body  was  measured 
using  a  cathetometer. 

The  five  pressure  ports  on  the  probe  were  connected  to  a  Scanivalve  Model  D 
pressure  scanner  capable  of  scanning  48  pressure  ports.  A  Hewlett  Packard  59306A 
Relay  Actuator  and  a  pressure  scanner  controller  actuated  the  pressure  scanner.  The 
pressure  from  each  port  of  the  pressure  scanner  was  read  by  a  Druck  Type  PDCR22 
differential  pressure  transducer  with  an  input  range  of  ±  1  psid. 

A  Pitot-static  tube  with  a  Datametric  Model  1173  manometer  measured  the  free- 
stream  dynamic  pressure.  A  Valdyne  Model  DB99  digital  barometer  measured  the  test 
section  static  pressure.  An  Instrulab  Model  1563  digital  thermometer  measured  the  free- 
stream  temperature.  The  signals  from  these  and  from  the  pressure  transducer  were 
scanned  by  a  Hewlett  Packard  3495  A  scanner  and  measured  by  a  Hewlett  Packard  3455 A 
digital  voltmeter.  A  Hewlett  Packard  9836  computer  controlled  the  data  acquisition 
system.  A  schematic  diagram  of  the  data  acquisition  and  control  system  is  shown  in 
Figure  5. 
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2.5  Hot-Wire  Anemometry  Instrumentation 
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Velocity  profiles  and  Reynolds  stress  components  in  the  boundary  layer  of  the 
inclined  prolate  spheroid  were  measured  using  a  DANTE C  55P61  x-wire  probe  in 
conjunction  with  two  Miller-type  (Miller,  1976)  integrated  circuit  constant  temperature 
hot  wire  anemometers  and  linearizers,  as  modified  by  Simpson,  et  al.  (1979)  The  x-wire 
probe  consists  of  two  platinum-plated  tungsten  wire  sensors  1 .25  mm  in  length  and  5  Mm 
in  diameter,  with  a  spacing  of  1  mm  between  the  sensor  wires.  The  geometry  of  the 
probe  is  shown  in  Figure  6.  An  overheat  ratio  of  1.7  was  used,  and  the  linearized 
calibrations  had  low'  dispersion  from  a  straight  line,  with  a  correlation  coefficient  greater 
than  0.999. 

To  obtain  all  six  Reynolds  stress  values  at  a  data  point,  an  x-wire  probe  must  be 
rotated  around  its  axis.  In  addition,  the  probe  should  be  aligned  within  10°  of  the  mean 
flow  direction  for  accuracy.  As  shown  in  Figure  7,  a  probe  support  with  rolling  and 
pitching  capabilities  was  designed  and  built  to  accommodate  these  requirements.  This 
support  was  manufactured  by  the  Virginia  Tech  Aerospace  and  Ocean  Engineering  Shop. 
When  mounted  on  the  Yj-Zt  traverse  described  previously  in  Section  2.3,  this  support 
allowed  the  x-wire  probe  to  be  pitched  from  0°  to  15°,  enabling  alignment  of  the  probe 
as  closely  as  possible  with  the  mean  velocity  for  a  given  profile.  A  Compumotor 
CXT32-39  stepper  motor  and  indexer  (Compumotor  Catalog,  1987)  controlled  by  an  IBM 
PS/2  Model  30  computer  through  an  RS-232  interface  was  utilized  to  roll  the  probe  along 
its  axis.  This  stepper  motor  allowed  highly  accurate  angular  displacement  of  the  probe. 
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with  25.000  microsteps  per  revolution  and  an  absolute  accuracy  of  0.0835°.  As  shown 
in  Figure  7,  a  guide  tube  with  two  struts  mounted  on  the  face  of  the  probe  support 
insured  that  the  probe  was  aligned  along  its  axis.  A  dial  indicator  and  a  cathetometer 
were  used  to  measure  the  probe  tip  displacement  from  centerline  to  adjust  the  alignment. 

As  showm  in  the  schematic  diagram  of  Figure  8,  the  voltage  signals  for  the  two 
wire  sensors  were  sampled  using  an  R.C.  Electronics  ISC- 16  A/D  board  mounted  in  an 
IBM  PS/2  Model  30  computer.  This  A/D  board  has  sixteen  channels  with  an  aggregate 
sampling  rate  of  up  to  1  Mhz  and  12-bit  accuracy  over  an  input  range  of  -10  to  +10 
volts.  A  Hewlett  Packard  3478A  multimeter  and  a  Tektronix  oscilloscope  were  used  to 
monitor  the  hot  wire  signals.  The  PS/2  Model  30  also  controlled  the  indexer  for  the 
stepper  motor  as  mentioned  previously.  The  Stability  Tunnel  Hewlett  Packard  computer 
system,  described  in  Section  2.4,  was  used  to  monitor  the  tunnel  status  and  obtain  the 
free-stream  conditions. 
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3.0  Experimental  Techniques  and  Data  Reduction 

3.1  Test  Conditions  and  Measurement  Locations 

The  investigations  were  performed  at  different  axial  locations  along  the  afterbody 
of  the  model.  These  locations  are  designated  by  X/L,  where  X  is  the  distance  measured 
from  the  nose  of  the  model  along  its  axis  of  symmetry,  and  L  is  the  length  of  the  model, 
1.37  meters  (54  inches).  Figure  9  shows  the  planes  where  the  measurements  were  made. 
These  are  X/L  =  0.6,  0.7,  0.8,  and  0.9.  At  each  axial  location,  data  were  obtained  in  a 
plane  parallel  to  the  YT-Z^  plane. 

For  the  outer  velocity  field  measurements,  the  model  was  inclined  at  angles  of 
attack  a  =  10Q  and  15°.  Measurements  were  made  at  X/L=0. 6,  0.7,  0.8,  and  0.9.  For 
all  the  locations,  free-stream  velocities  of  approximately  15.2  mps  (50  fps)  and  45.7  mps 
(150  fps)  were  used,  giving  Reynolds  numbers  based  on  the  length  of  the  model  of 
1.3x  106  and  4.0 x  106,  respectively,  at  a  temperature  of  25°C. 
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Boundary  layer  mean  velocity  and  turbulence  measurements  were  obtained  for  a 


Reynolds  number  of  4. Ox  106  at  a=  10°  and  X/L  =  0.7  and  0.8. 


3.2  Mean  Velocity  Field  Measurements 

Outer  mean  velocity  field  measurements  in  the  YT-Z ^  plane  at  each  axial  location 
of  the  inclined  prolate  spheroid  were  made  using  the  yawhead  probe  described  in  Section 
2.4.  At  each  data  point,  the  five  ports  were  scanned  and  the  pressure  from  each  port  was 
averaged  over  25  samples.  A  computer  program  used  by  Choi  and  Simpson  (1987) 
controlled  the  scanning  operation.  This  program  also  performed  online  data  reduction 
using  the  yawhead  probe  calibration  and  data  reduction  method  of  Sung  (1985).  This 
method  is  explained  in  Appendix  A.  The  free-stream  temperature,  dynamic  pressure,  and 
barometric  pressure  were  recorded  for  each  data  point.  The  velocity  components  were 
first  found  in  probe  coordinates  and  then  transformed  into  tunnel  coordinates,  taking  into 
account  probe  misalignment  with  the  XT-axis.  Probe  misalignment  angles  are  showm  in 
Figure  10  along  with  their  relation  to  tunnel  coordinates. 

The  sensitivity  of  the  differential  pressure  transducer  was  found  using  a  three  point 
calibration.  An  AMTEK  Model  MK  100  pneumatic  pressure  tester  was  used  to  apply 
pressures  of  0  psid  and  ±  1  psid  to  the  transducer,  and  the  transducer  voltage  output  was 
assumed  linear  in  this  range.  During  the  data  acquisition,  a  sixth  port  on  the  Scanivalve 
pressure  scanner  was  fed  1  psi  of  pressure  to  monitor  the  calibration  of  the  transducer. 
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3.3  Velocity  Profiles  and  Turbulence  Measurements 


Velocity  profiles  and  Reynolds  stress  quantities  near  the  surface  of  the  inclined 
prolate  spheroid  were  obtained  in  the  Yt-Zt-  plane  at  two  axial  locations  using  the  x-wire 
probe  supported  by  the  rolling  and  pitching  mount  described  in  Section  2.5.  Figure  6 
shows  the  probe  coordinate  system  and  the  definitions  of  roll  angle  4r  and  wire  angles  0j 
and  0:.  To  obtain  all  six  Reynolds  stress  quantities  and  all  three  mean  velocity 
components,  the  probe  was  rolled  to  the  following  positions: 

(a)  U-V  plane  (ijr  =  0°) 

(b)  +45°  plane  (ijr  =  45°) 

(c)  U-W  plane  (t|/  =  90°) 

(d)  +135°  plane  (\|r  =  135°) 

(e)  U-V  plane  (tjr  =  180°) 

A  description  of  the  data  reduction  method,  equations  used,  and  derivations  is  given  in 
Appendix  B. 

Initial  calibrations  of  the  x-wire  probe  were  made  using  a  TSI  Model  1 125 
calibrator  at  various  probe  orientations  with  respect  to  the  flow  to  obtain  tangential  and 
normal  cooling  factors.  A  description  of  the  choice  for  these  factors  is  given  in  Appendix 
B.  Before  and  after  each  run.  the  probe  was  placed  in  the  tunnel  free  stream  and  aligned 
in  the  Xx  direction,  and  in-tunnel  calibrations  were  made  at  various  flow  speeds.  As 


Experimental  Techniques  and  Data  Reduction 


13 


mentioned  previously,  each  linearized  calibration  had  a  correlation  coefficient  greater  than 
0.999. 


Before  acquiring  data,  the  probe  was  pitched  to  the  desired  angle,  based  on  the 
previous  mean  flow  measurements,  to  align  the  probe  as  closely  as  possible  to  the  mean 
velocity  for  a  given  profile.  The  probe  was  adjusted  to  the  proper  axial  (Xx)  location 
along  the  model.  Setting  the  pitch  angle  of  the  probe  and  positioning  the  probe  at  the 
proper  Xx  location  were  accomplished  using  the  probe  support  with  adjustable  pitch  angle 
capabilities  described  in  Section  2.5  and  shown  in  Figure  7.  A  caihetometer  was  used 
to  measure  the  pitch  angle  of  the  probe  with  respect  to  the  Xx  axis.  The  Yx-Z^  traverse 
was  used  to  place  the  probe  in  the  correct  position  for  a  given  profile.  During  the 
measurements,  the  probe  was  moved  to  each  data  point  in  the  profile  using  the  Yr-2^ 
traverse.  A  cathetometer  was  used  to  measure  the  location  of  the  probe  relative  to  the 
model.  At  each  data  point,  the  probe  was  rolled  about  its  axis  to  the  angles  described 
above.  For  each  angle,  the  voltage  pairs  were  sampled  at  820  Hz  to  obtain  16384 
samples  over  a  twenty  second  period.  Mean  and  RMS  quantities  of  the  voltages  were 
stored  on  disk.  In  addition,  the  tunnel  free-stream  dynamic  pressure,  temperature,  and 
barometric  pressure  were  recorded  for  each  data  point. 

The  data  were  reduced  using  the  method  outlined  in  Appendix  B  to  obtain  mean 
velocity  components  and  Reynolds  stresses  at  each  point  in  probe  coordinates,  and  a 
tensor  transformation  was  used  to  find  these  values  in  tunnel  coordinates.  The  method 
of  Larsen  and  Busch  (1974),  described  in  Appendix  B,  was  used  to  correct  for 
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temperature  differences  between  the  calibrations  and  the  acquired  data.  The  computer 
program  listed  in  Appendix  B  was  used  to  perform  the  data  reduction. 

3.4  Uncertainty  Analysis 

Uncertainties  in  the  yawhead  probe  measurements  include  uncertainties  in 
measuring  the  probe  misalignment  with  respect  to  the  XT-axis.  Other  uncertainties 
include  those  in  velocity  vector  field  measurements  made  using  the  yawhead  probe.  The 
following  RMS  errors  are  given  by  Sung  (1985)  for  this  probe: 


Pitch  angle  a 

±  3.93° 

Yaw  angle  6 

±  2.36° 

Stagnation  pressure  P, 

±  4.82% 

Dynamic  pressure  q 

±  5.33% 

Total  Velocity  Ul0, 

±  2.73% 

The  effects  of  these  uncertainties  on  the  measured  mean 

were  estimated,  and  the  following  values  were  found: 

velocities  in  tunnel  coordinates 

Mean  Velocity  U 

±  0.035U 

Mean  Velocity  V 

±  0.045U 

Mean  Velocity  W 

±  0.050U 

Experimental  Techniques  and  Data  Reduction 

15 

Appendix  B  briefly  details  the  agreement  among  the  estimates  for  each  of  the 
quantities  obtained  using  the  x-wire  probe.  In  addition,  various  sources  of  uncertainty 
exist  that  must  be  considered,  including  the  following: 

1 .  Probe  position 

2.  Probe  pitch  angle  misalignment 

3.  Angle  of  attack  of  model 

4.  Hot-wire  instrumentation  and  calibration 

5.  Resolution  of  voltage  signals 

6.  Probe  geometry 

7.  Tunnel  free-stream  temperature 

8.  Tunne.  flee-stream  dynamic  pressure 

9.  Tunnel  static  pressure 

An  attempt  was  made  to  determine  the  magnitude  of  the  uncertainty  associated  with  each 
source  of  error.  An  additional  source  of  error  exists  due  to  the  1  mm  spacing  between 
the  two  sensors  on  the  x-wire  probe.  Because  the  probe  is  rolled  about  its  longitudinal 
axis  to  obtain  all  six  Reynolds  stress  components,  the  data  are  acquired  for  a 
measurement  volume,  instead  of  a  measurement  point.  Therefore,  some  uncertainties  in 
the  measured  quantities  exist  due  to  this  uncertainty  in  the  actual  measurement  location. 

The  overall  effects  of  the  uncertainties  listed  above  on  the  mean  velocities  and  the 
turbulence  quantities  were  investigated  using  a  perturbation  computer  program 
incorporating  the  methods  of  Kline  and  McClintock  (1953)  and  Moffat  (1982).  The 
following  estimates  for  the  uncertainties  m  the  measured  quantities  were  found: 
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• 

Mean  Velocity  U 

±  0.033U 

• 

Mean  Velocity  V 

±  0.035U 

Mean  Velocity  W 

±  0.040U 

• 

Reynolds  Stress  u 2 

+  0.090  u2 

—  max 

Reynolds  Stress  v2 

*  0.090  tT2  T 

—  max 

• 

Reynolds  Stress  H': 

*  0.12^ 

• 

Reynolds  Stress  uv 

+  0.  16uv 

—  max 

Reynolds  Stress  uw 

+  0.17uv 

—  max 

• 

Reynolds  Stress  tw 

±  0.27  uv 

max 

For  the  uncertainties  in  the  Reynolds  stress  values. 

u2  and  uv 

max  max 

refer  to  the 

• 

maximum  values  in  each  measurement  profile. 

• 

• 

• 
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• 

4.0  Experimental  Results  and  Discussion 

4.1  Mean  Velocity  Field  Measurements 

Mean  velocity  field  measurements  were  made  in  vertical  planes  along  the  leeward 
side  of  the  body  at  X/L=0.6.  0.7.  0.8.  and  0.9,  measured  along  the  centerline  axis  of 
the  model.  The  model  was  placed  at  angles  of  attack  a  of  10°  and  15°  with  respect  to 
the  wind  tunnel  free  stream,  and  data  were  obtained  for  Reynolds  numbers  of  1.3  x  106 
and  4.0  x  106.  The  measurements  were  made  with  a  grid  spacing  of  1.27  cm  (0.5  inches) 
in  the  YT  and  Zj  directions  in  each  plane,  from  Yx  =  -2.54  cm  to  21.57  cm  and  from 
Zt  =  -15.23  cm  to  2.54  cm.  Components  U,  V,  and  W  of  the  velocity  vector  at  each  data 
point  were  found  in  tunnel  coordinates. 

Before  discussing  these  results,  it  is  worthwhile  to  consider  the  surface  oil-flow 
visualization  experiments  conducted  by  Ahn  (1990)  on  the  same  6:1  prolate  spheroid 
model  at  various  conditions.  The  results  of  these  oil-flow  experiments  reveal  skin  friction 
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lines  and  regions  of  converging  streamlines  on  the  surface,  and  these  converging 
streamlines  can  indicate  the  locations  of  surface  flow  separation. 

These  surface  oil-flow  experiments  cannot  conclusively  show  whether  the  flow 
over  the  6:1  prolate  spheroid  at  Re=1.3xl06  and  a  =  10°  and  15°  is  laminar  or 
turbulent.  How-ever,  the  results  of  surface  skin  friction  measurements  on  a  6:1  prolate 
spheroid  performed  by  Kreplin.  ei  al.  (1980)  indicate  that  at  Re=1.6xl06  and  a=10° 
the  flow'  is  laminar  up  to  separation.  Thus,  in  the  present  study  at  Re  =  1 . 3  x  1 06  and 
a=  10°  and  15°.  the  flow  over  the  body  is  possibly  laminar  up  to  separation. 

Results  of  surface  oil-flow  visualizations  for  Re=4.0x  106  and  a  =  10°  are  shown 
in  Figures  1 1  and  12.  For  these  conditions,  transition  occurs  along  the  forebody  and  the 
flow  becomes  turbulent.  Surface  skin  friction  lines  extending  from  the  windward  side  of 
the  forebody  converge  towards  skin  friction  lines  from  the  leeward  side.  At  X/L=0.6, 
this  convergence  occurs  between  <p  =  35°  and  40°;  at  X/L=0.7,  it  occurs  between  tp  =  25  ° 
and  30°;  at  X/L=0.8,  it  occurs  between  <p=l5°  and  20°;  and  at  X/L=0.9,  the 
convergence  occurs  below  <p=10°.  Between  X/L=0.8  and  0.9,  primary  separation 
begins  along  the  line  of  converging  skin  friction  lines.  Aft  of  X/L=0.9,  secondary 
separation  is  evident. 

For  a=  15°  and  Re=4.0x  106,  oil-flow  results  are  shown  in  Figures  13  and  14. 
For  these  conditions,  transition  occurs  around  X/L  =  0.3  and  the  flow  becomes  turbulent. 
Regions  of  converging  skin  friction  lines  are  evident  at  X/L  =  0.6  between  <p=40°  and 
45°,  at  X/L  =  0.7  between  <p  =  30°  and  35°,  at  X/L  =  0.8  between  ip  =  20°  and  25°.  and 
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at  X/L=0.9  around  <p  =  0°.  Primary  separation  occurs  farther  forward  along  the  body 
than  for  a=10°,  starting  at  X/L  =  0.5. 

Based  upon  discussions  of  the  topology  of  separated  flow,  such  as  those  by 
Lighthill  (1963)  and  Peake  and  Tobak  (1982),  the  separating  flow  along  this  body  can  be 
considered  as  a  sheet  of  fluid  that  comes  out  from  the  surface  in  the  regions  of 
converging  skin  friction  lines  indicated  by  the  oil-flow  experiments.  Along  the  surface 
of  the  body,  this  sheet  rolls  up  by  turning  in  toward  the  leeward  side  of  the  body,  down 
toward  the  surface,  then  away  from  the  leeward  side,  as  indicated  in  Figure  15.  In  the 
context  of  this  discussion,  the  phrase  "vortical"  refers  to  this  roll-up  of  the  separating 
sheet  of  fluid. 

When  viewed  along  the  longitudinal  axis  of  the  body,  vector  plots  of  the 
secondary-  velocity  components  can  reveal  the  nature  of  the  flow  field  as  the  separating 
fluid  turns  and  rolls  over  in  the  manner  previously  described.  Therefore,  the  data  were 
transformed  by  a  rotation  about  the  Zy  axis  of  -a  to  a  body-fixed  coordinate  system  to 
obtain,  in  effect,  a  view  along  the  longitudinal  body  axis.  This  body-fixed  system  is 
denoted  by  the  axes  XB,  YB,  and  ZB;  with  XB  corresponding  to  the  longitudinal  axis  of 
the  body.  The  rollover  of  the  flow  becomes  more  apparent  in  the  resulting  secondary 
velocity  vector  plots,  as  can  be  seen  by  comparing  Figure  16  to  Figure  32.  For 
Re  =  1.3x  106,  Figures  17  to  20  show  the  results  for  a=  10°  and  Figures  21  to  24  show 
the  results  for  a=  15°.  For  Re  =  4.0x  1 06,  Figures  25  to  28  show  the  results  for  a=  10° 
and  Figures  29  to  32  show  the  results  for  a  =  15°. 
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For  a  =  10°  and  Re  =  1 .3  x  106.  the  secondary  velocity  vectors  show  some  turning 
of  the  flow,  as  indicated  by  the  change  in  the  direction  of  the  V  component  of  velocity. 
The  magnitude  of  turning  increases  and  spreads  further  from  the  surface  as  the  flow 
progresses  toward  the  tail  of  the  body. 

For  a=  10°  and  Re=4.0x  106.  the  turning  of  the  flow  is  again  shown  as  described 
above.  However,  the  magnitude  of  the  turning  increases  in  such  a  manner  that,  starting 
at  X/L  =  0.8  and  continuing  at  X/L  =0.9.  a  region  of  definite  vortical  flow  appears, 
although  the  location  of  its  center  is  difficult  to  ascertain. 

For  Re=1.3xl06,  the  turning  of  the  flow  at  each  axial  location  occurs  further 
from  the  surface  at  a  =  15°  than  at  a=  10°.  Vortical  flow  encompassing  a  large  region 
appears  starting  at  X/L=0.8.  The  small  magnitudes  of  the  V  and  W  components  indicate 
the  weak  nature  of  this  flow. 

For  a=15°  and  Re  =  4.0xl06,  a  distinct  vortical  pattern  appears  close  to  the 
surface  at  X/L  =  0.6.  and  the  magnitudes  of  the  secondary  velocities  increase  and  the 
region  of  vortical  flow  grows  as  the  flow  progresses  aft  along  the  body  to  X/L=0.9. 
This  vortical  flow  appears  to  be  strong,  considering  the  magnitude  of  the  secondary 
velocities.  Again,  the  location  of  the  center  of  the  region  is  difficult  to  judge,  although 
the  center  appears  to  be  moving  away  from  the  surface  at  X/L =0.9. 

More  detailed  information  about  the  development  of  the  vortical  flow  on  the 
leeward  side  of  the  body  for  each  test  condition  can  be  seen  by  examining  contours  of 
streamwise  vorticity  at  each  axial  location.  In  the  presence  of  a  vortex,  the  location  and 
size  of  the  vortex  can  be  found  by  studying  any  concentric  contours  that  might  appear. 
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and  the  core  of  the  vortex  can  be  located  by  finding  the  position  of  the  maximum 
vorticity  within  such  concentric  contours  (Pauley  and  Eaton,  1988).  In  addition,  a  region 
or  layer  of  strong  streamwise  vorticity  growth  can  be  seen  in  closely  spaced  contour  lines 
across  which  occurs  a  large  increase  in  vorticity.  Such  a  growth  layer  coincides  with  a 
line  of  separation  in  the  flow  field,  as  described  by  Barberis  and  Chanetz  (1986).  These 
curves  of  constant  vorticity  can  give  a  good  image  of  a  sheet  of  fluid  that  separates  from 
the  body  and  winds  up  to  form  vortical  flow.  By  studying  vorticity  contours  for  the 
present  conditions,  some  conclusions  can  be  made  concerning  the  size,  location,  and 
strength  of  any  vortical  flow  regions. 

Measurements  of  the  secondary  velocity  components  V  and  W  in  tunnel 
coordinates  allow,  at  each  data  point,  the  determination  of  the  vorticity  in  the  Xx 
direction,  given  by: 


dV' 

dZ, 


The  local  gradients  in  the  YT  and  Z \  directions  were  found  by  fitting  a  parabola  through 
five  consecutive  data  points  and  calculating  the  derivative  at  the  central  point.  A  possible 
source  of  error  in  the  values  calculated  for  the  gradients  exists  due  to  the  spacing  of  the 
data  points.  The  1.27  cm  spacing  may  be  too  large  to  accurately  capture  the  gradients 
in  each  direction.  No  attempt  was  made  to  estimate  the  possible  error.  The  resulting 
contours  of  nx/U„  are  shown  in  Figures  33  to  48. 

The  results  for  a  =  10°  and  Re=1.3xl06  are  shown  in  Figures  33  to  36.  At 
X/L=0.6,  the  only  noticeable  feature  is  a  region  between  tp  =  35°  to  ip  =  50°  close  to  the 
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surface  in  which  the  vorticity  nx/U„  increases  from  0.0/cm  to  0. 10/cm.  Note  that  tp  is 
the  angle  measured  from  the  -Zj  axis  as  shown  in  Figure  4.  At  X/L  =  0.7,  this  same  type 
of  region  is  seen.  In  addition,  between  <p  =  45°  and  <p  =  75°.  the  contours  spread  out  from 
the  surface,  indicating  diffusion  of  vorticity  into  the  flow  field  surrounding  the  body. 
This  spreading  is  even  more  evident  at  X/L  =  0.8  between  <p=40°  and  <p  =  80°.  The 
vorticity  nx/U„  increases  from  0.0/cm  to  a  maximum  of  0.10/cm  between  <#>=45°  and 
<P  =  65°.  At  X/L=0.9,  contours  are  almost  concentric  at  higher  ip's,  increasing  from 
0.0/cm  at  the  edges  to  0.08/cm  towards  the  center.  However,  the  maximum  value  occurs 
at  approximately  <p  =  35°.  in  a  region  of  vorticity  growth  from  0.0/cm  to  0.14/cm  close 
to  the  surface. 

These  vorticity  contours  for  a=10°  and  Re=1.3xl06  show-  the  diffusion  of 
vorticity  out  into  the  flow  field  as  indicated  by  regions  of  vorticity  that  originate  close 
to  the  body  and  spread  out  away  from  the  surface  as  the  flow  progresses  aft  along  the 
body.  The  maximum  vorticity  nx/U„  is  fairly  constant  at  0. 10/cm  until  X/L=0.9,  when 
it  reaches  0.14/cm,  indicating  a  slight  increase  due  to  the  early  stages  of  the  rollup  of 
separating  fluid  at  the  tail  of  the  body. 

Figures  37  to  40  show  the  nx/U„.  vorticity  contours  for  a=  10°  and  Re=4.0x  106. 
At  X/L  =  0.6,  the  vorticity  nx/U_  goes  from  0.0/cm  to  0.06/cm  at  <p  =  70°  to  75°.  This 
is  at  a  higher  ip  than  the  low'er  Re  case  at  the  same  axial  location  and  angle  of  attack. 
At  X/L  =  0.7,  the  vorticity  increases  from  0.0/cm  to  0. 14/cm  in  a  region  between  <p=40° 
and  <p  =  50°.  The  contours  are  spreading  out  away  from  the  surface  between  <p=60°  and 
tp  =  85°.  At  X/L=0.8,  the  vorticity  nx/U„  increases  from  0.0/cm  to  a  maximum  of 
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0.16/cm  across  a  layer  in  the  same  <p  locations  as  at  X/L=0.7.  Again,  the  contours  are 
spreading  out  away  from  the  body  between  tp  =  60°  and  <p  =  85°.  Contours  that  are  almost 
concentric  are  evident  at  X/L=0.9,  indicating  vortical  flow  centered  between  tp  =  55°  and 
ip  =  60°  and  approximately  2  cm  out  from  the  surface.  The  maximum  vorticity  nx/U„ 
is  greater  than  0. 16/cm  at  the  center  of  this  region. 

These  contours  show  that  the  separation  results  in  vortical  flow'  centered  close  to 
the  body  at  X/L=0.9.  The  contours  appear  to  take  on  an  elliptical  shape.  The  maximum 
vorticity  nx/U„at  each  axial  location  increases  from  0.06/cm  at  X/L=0.6  to  greater  than 
0. 16/cm  at  X/L=0.9. 

For  a  =15°  and  Re=1.3xl06.  the  constant  vorticity  contours  are  shown  in 
Figures  41  to  44.  At  X/L=0.6,  the  values  increase  from  0.0/cm  to  0.16/cm  between 
ip  =  35°  and  <p  =  50°.  At  higher  ip’s.  the  contour  lines  spread  out  from  the  body.  The 
vorticity  nx/U„  increases  from  0.0/cm  to  0.20/cm  between  tp-5°  and  45°  at  X/L=0.7. 
Again,  spreading  of  the  contours  is  seen  at  higher  ip's.  At  X/L=0.8  the  vorticity  nx/U_ 
increases  from  0.0/cm  to  0.20/cm  in  a  region  centered  at  <p=40°.  At  X/L=0.9.  the 
values  increase  from  0.0/cm  to  0.32/cm  between  <p's  of  0°  and  35°.  An  interesting 
feature  to  note  is  that  contour  lines  of  0.04/cm  and  0.08/cm  have  wrapped  around  and 
moved  close  to  the  surface  between  <p=45°  and  <p  = 90°.  This  was  not  seen  in  the  case 
with  a=  10°  and  Re=  1 .3  x  106.  These  wrapped  contours  indicate  that  the  separation  fluid 
sheet  has  started  to  roll  up,  resulting  in  vortical  flow  centered  between  <p  =  60°  and  75°. 
with  a  maximum  vorticity  in  the  center  of  approximately  0.10/cm.  The  spread  of  the 
contour  lines  indicate  diffusion  of  vorticity  over  a  large  area  away  from  the  body. 
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As  these  results  for  the  case  of  a=15°  and  Re=1.3xl06  show,  the  maximum 
vorticity  nx/U„  increases  from  0.16/cm  to  0.32/cm  between  X/L=0.6  and  X/L  =  0.9, 
and  in  all  four  axial  locations  this  maximum  vorticity  occurs  close  to  the  surface  at  <p's 
less  than  45°.  Contour  lines  indicate  diffusion  of  streamwise  vorticity  over  a  large  region 
away  from  the  body  at  <p’s  greater  than  45°.  Weak,  rapidly  diffusing  vortical  flow  due 
to  the  rollup  of  the  separating  fluid  appears  to  have  formed  by  X/L=0.9. 

The  results  for  the  case  of  a  =  15°  and  Re  =  4.0xl06  are  shown  in  Figures  45  to 
48.  At  X/L=0.6,  a  maximum  vorticity  nx/U„  of  0.16/cm  is  seen  at  approximately 
ip  =  55°.  Between  <p  =  60°  and  <p  =  85°,  contour  lines  spread  out  away  from  the  surface. 
At  X/L  =  0.7 .  a  maximum  vorticity  of  0.25/cm  occurs  in  a  contour  that  is  almost  closed, 
indicating  early  stages  of  the  formation  of  vortical  flow  centered  at  approximately 
<p  =  55°.  In  the  region  above  ip  =  60°,  the  contours  are  spread  out  slightly  compared  to 
those  at  X/L=0.6.  Definite  concentric  contours  are  seen  at  X/L=0.8,  indicating  vortical 
flow  centered  2  cm  from  the  surface  at  ip  =  55°  to  65°.  The  maximum  vorticity  nx/U„ 
at  the  center  of  these  contours  is  0.30/cm.  Note  that  the  vorticity  contours  have  not 
spread  out  away  from  the  body  any  more  than  at  X/L  =  0.7,  indicating  that  the  vortical 
flow  contains  the  vorticity  and  slows  its  diffusion.  At  X/L  =  0.9,  the  center  has  moved 
away  from  the  body,  as  was  also  shown  previously  by  the  secondary  velocity  vectors. 
The  vortical  flow  is  centered  3  cm  from  the  surface  at  <p  =  60°  and  the  contours  have  a 
definite  elliptical  shape.  The  strength  of  the  vortical  flow  is  indicated  by  the  maximum 
vorticity  nx/U„of  0.36/cm  at  its  center.  Again,  diffusion  of  vorticity  appears  to  be 
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prevented  in  the  presence  of  this  flow,  as  seen  by  the  fact  that  contour  lines  are  contained 
without  spreading. 

As  these  results  show  for  a=15°  and  Re=4.0xl06,  strong  vortical  flow  has 
formed,  originating  between  X/L=0.6  and  X/L  =  0.7  and  gaining  in  strength  to  X/L  =  0.9 
as  indicated  by  the  increase  in  maximum  vorticity  nx/LL  from  0.16/cm  to  0.30/cm. 

Overall,  the  results  from  the  secondary  velocity  vectors  and  the  contours  of 
vorticity  in  the  XT  direction  indicate  that  for  Re=  1.3x  106  at  both  angles  of  attack,  any 
vortical  formation  is  very  weak  and  rapidly  diffusing.  In  contrast,  both  cases  at 
Re  =  4.0x  10*  show  the  formation  of  vortical  flow  of  increasing  strength. 


4.2  Mean  Velocity  Profiles  and  Turbulence  Measurements 


Mean  velocity  profiles  and  turbulence  measurements  were  made  at  X/L  =  0.7  and 
X/L=0.8  at  a  =  10°  and  Re  =  4.0xl06.  Data  were  obtained  in  vertical  and  horizontal 
traverses  in  the  YT  and  Zj-  directions  in  each  plane  as  shown  in  Figures  49  and  50. 

In  Figures  51  to  57,  mean  velocity  components  obtained  using  the  x-wire  probe 
are  compared  to  mean  velocity  components  obtained  using  the  five-hole  yawhead  probe 
in  vertical  profiles  at  the  same  locations.  As  these  comparisons  show,  the  values  of 
U/Ure  agree  to  within  approximately  4%.  In  most  cases,  the  values  of  V/Ure  and  W/Ure 
agree  to  within  10%.  Consideration  of  these  comparisons  appears  to  show  that  the 
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uncertainties  in  the  mean  velocities  discussed  in  Section  3.4  are  reasonable,  and  any 
disagreements  may  possibly  be  attributed  to  interference  from  the  probes  and  supports. 

To  discuss  the  results  and  compare  with  work  performed  by  other  researchers,  the 
data  were  linearly  interpolated  to  obtain  profiles  extending  radially  outward  from  the 
surface  of  the  body.  As  shown  in  Figures  49  and  50,  these  profiles  are  located  from 
ip  =  20°  to  <p  =  90°  at  each  axial  location. 


4.2.1  Velocity  Profile  Measurements 


Mean  velocity  profile  measurements  are  presented  in  a  local  body  coordinate 
system  shown  in  Figure  58.  The  origin  of  this  system  is  on  the  surface  of  the  model, 
with  the  Xb  axis  aligned  with  the  axis  of  the  model,  the  Yb  axis  extending  radially 
outward  from  the  surface,  and  the  7^  axis  tangential  to  the  surface.  The  data  were 
transformed  by  rotating  around  the  Zj  axis  by  -a  and  then  rotating  around  the  X  axis  by 
-(90°-<p)  to  the  Xb-Yb-Zi,  coordinate  system  for  each  profile.  In  this  coordinate  system, 
for  each  data  point,  Yb  is  the  distance  from  the  model  surface  and  Zb=0.0.  This 
coordinate  system  was  chosen  to  make  comparisons  with  the  data  of  Meier,  el  al.  (1984a), 
because  in  this  system  W  is  the  spanwise  (circumferential)  velocity  component  and  a 
change  in  the  sign  of  W  corresponds  to  a  spanwise  flow  reversal.  The  velocities  were 
nondimensionalized  by  the  local  free  stream  velocity,  Ure.  This  local  free  stream  velocity 
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for  each  radial  profile  was  determined  by  plotting  the  resultant  velocity 
Ur=(U:  +  V:  +  W2)1'*  versus  Yb  and  determining  Ur  at  the  edge  of  the  boundary  layer. 
This  was  considered  to  be  Urc  for  that  profile. 

The  resulting  velocity  profiles  are  shown  in  Figures  59  to  74.  Due  to  the  size  of 
the  x-wire  probe  used  to  obtain  the  data,  measurements  could  not  be  made  closer  than 
approximately  0.25  cm  from  the  surface.  Therefore,  almost  all  of  the  data  are  in  the 
outer  region  of  the  boundary  layer  and  no  measurements  were  made  in  the  inner  region 
or  laminar  sublayer.  This  is  reflected  in  the  velocity  profiles. 

The  boundary  layer  thickness  at  each  radial  location  was  estimated  by  calculating 
<5W,  given  by  Yh  where  Ur=0.99Ure.  The  displacement  thickness  6*  and  the  momentum 
thickness  0.  defined  by: 


can  give  an  indication  of  the  growth  of  the  boundary  layer  over  the  leeward  side  of  the 
body.  To  estimate  these  values,  the  velocity  profiles  were  extended  to  the  surface  of  the 
body  using  the  following  log  law  for  a  turbulent  flat-plate  boundary  layer  taken  from 
Schetz  (1984): 
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The  three  data  points  closest  to  the  surface  in  each  profile  were  used  to  determine  the 
values  of  A.  C,  and  u.  for  that  profile.  The  displacement  thickness  and  momentum 
thickness  were  then  calculated  using  the  trapezoidal  rule  to  evaluate  the  integrals.  The 
shape  factor  for  each  profile.  H  =  <5*/0.  was  also  calculated.  Results  are  given  in  Table 
2  and  Table  3.  Although  this  log  lav.  is  not  exactly  valid  for  the  present  flow,  the 
method  of  extending  the  profiles  to  the  surface  provides  better  estimates  of  the 
displacement  thickness  and  momentum  thickness  than  those  obtained  when  neglecting 
such  extensions. 

Figure  75  shows  the  circumferential  development  of  the  boundary'  layer  given  by- 
displacement  thickness  6*  versus  tp  for  X/L— 0.7  and  X/L=0.8.  This  figure  also  shows 
the  results  from  Meier,  et  al.  (1984a)  for  Re  =  7.2xl06  at  X/L=0.73  and  X/L=0.64. 
The  effects  of  Reynolds  number  on  the  growth  of  the  boundary  layer  are  clearly  evident. 
At  the  higher  Reynolds  number,  the  displacement  thickness  is  greater  and  the  location  of 
the  maximum  displacement  thickness  is  shifted  towards  the  leeward  side  of  the  body.  A 
possible  explanation  for  the  larger  displacement  thickness  at  the  higher  Reynolds  number 
is  given  by  considering  the  locations  of  transition  in  each  case.  At  Re  =  7.2xl06, 
boundary  layer  transition  begins  to  occur  in  an  axial  location  closer  to  the  nose  of  the 
body  than  at  Re=4.0xl0f\  For  the  higher  Reynolds  number,  this  earlier  transition 
allows  the  turbulent  boundary  layer  to  grow-  more  before  reaching  the  locations  of 
separation. 
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The  profiles  of  circumferential  crossflow  W/Ure,  in  Figures  59  to  74.  indicate  a 
change  of  sign  of  the  crossflow  close  to  the  surface  of  the  body  at  <p  =  50°  for  X/L  =  0.7 
and  at  slightly  over  ip=40°  for  X 1  =0.8.  taking  into  account  that  the  measurements  were 
not  made  closer  than  0.25  cm  from  the  surface.  Data  from  Meier,  ei  al.  (1984a)  show 
that,  at  the  higher  Reynolds  number,  this  negative  crossflow  occurs  at  <p  =  35°  for 
X’L  =  0.73  and  at  <p  =  50°  for  X/L  =  0.64.  Thus,  the  higher  Reynolds  number  causes 
negative  crossflow  to  occur  further  away  from  the  leeward  side  at  a  given  axial  location. 
This  negative  crossflow  coincides  roughly  with  a  substantial  increase  in  the  displacement 
thickness  at  the  locations  that  it  occurs,  as  shown  in  Figure  75.  Meier,  ei  al.  (1984a) 
report  that  the  movement  of  the  location  of  negative  crossflow  awray  from  the  leeward 
side  is  accompanied  by  a  shift  in  the  location  of  the  maximum  displacement  thickness. 
The  present  data  appear  to  support  this.  The  change  in  sign  of  the  crossflow  W/Ure  near 
the  surface  does  not  indicate  the  location  of  separation,  because  the  location  of  negative 
crossflow  depends  upon  the  choice  of  coordinate  system  (Meier,  et  al.,  1983b). 

Although  at  <p  =  90c  the  crossflow  velocity  W/Ure  should  be  approximately  zero, 
the  data  at  X/L  =  0.7  in  Figure  66  give  non-zero  values.  However,  most  of  these  values 
are  within  the  uncertainty  range  for  W  given  in  Section  3.4 

In  addition  to  showing  the  velocity  profiles.  Figures  59  to  74  present  hodograph 
plots  that  show  the  variation  of  tne  direction  of  U-W  crossflow  velocity  vectors  between 
the  edge  of  the  boundary  layer  and  the  surface,  in  the  Z ^  direction.  This  variation,  or 
skewing,  has  significance  because  in  a  three-dimensional  boundary  layer,  fluid  can  move 
laterally,  as  opposed  to  a  two-dimensional  layer  in  which  the  fluid  is  constrained  on  time 
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average  to  a  plane.  Investigating  such  skewing  is  difficult  in  the  presence  of  large 
rotation  of  the  flow  because  the  distinction  between  the  outer  inviscid  fluid  and  the 
boundary  layer  is  blurred  (Barberis  and  Chanetz,  1986).  However,  the  present  data  were 
obtained  in  a  region  without  large  fluid  rotation,  and  studying  this  variation  in  the 
direction  of  U-W  crossflow  vectors  can  be  meaningful.  The  velocity  vector  at  the  surface 
is  in  a  direction  parallel  to  the  surface  stress,  and  the  variations  show  how  the  direction 
of  the  surface  stress  vectors  lag  when  compared  with  the  direction  of  the  local  free-stream 
velocity  (Townsend,  1976). 

Figure  76  and  77  show  profiles  of  the  flow  angle,  B  =  arctan(W/U).  in  local  body 
coordinates,  for  all  radial  profiles  at  X/L=0.7  and  0.8.  These  give  an  indication  of  the 
variation  of  the  direction,  or  skewing,  of  the  flow  through  the  layer.  The  amount  of 
skewing  in  each  radial  profile  can  be  determined  by  finding  the  difference  between  the 
flow  angle  of  the  local  free  stream,  Bfs,  and  the  flow  angle  of  the  velocity'  at  the  surface, 
Bsurf.  This  difference  is  given  by  AB  =  B(s-Bsurf.  However,  because  measurements  were 
not  obtained  close  to  the  surface  in  the  laminar  sublayer,  the  present  data  can  give  only 
an  estimate  of  AB.  The  results  are  showm  in  Figure  78,  along  with  data  from  Meier,  et 
a/. (1984a)  taken  at  Re  =  7.2x  106.  A  higher  Reynolds  number  appears  to  cause  greater 
skewing  through  the  boundary  layer,  in  effect  causing  greater  lag  in  the  surface  stress. 
The  maximum  skewing  for  all  the  axial  locations  occurs  at  <p  =  60°;  thus,  this  appears  to 
be  independent  of  the  location  of  negative  crossflow  W/Urc  and  maximum  displacement 
thickness.  A  comparison  between  the  location  of  maximum  skewing  and  the  location  of 
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maximum  <5*  for  a  given  axial  location  shows  that  the  maximum  skewing  occurs  closer 
to  the  leeward  side. 


4.2.2  Turbulence  Measurements 

Turbulence  measurements  are  presented  in  a  local  free-stream  velocity  coordinate 
system  shown  in  Figure  79.  The  origin  of  this  system  is  on  the  surface  of  the  model, 
with  the  Xfs  axis  aligned  with  the  local  free-stream  velocity  Ure  at  the  edge  of  the 
boundary  layer,  the  Yfs  axis  extending  radially  outward  from  the  surface,  and  the  Z,s  axis 
completing  the  right-hand  rule.  With  the  components  of  Urc  denoted  as  Ue.  Ve,  and  We. 
the  data  were  transformed  from  tunnel  coordinates  by  rotating  around  the  Z^  axis  by 
afs  =  tan'1(Vc/Uc).  rotating  around  the  Y  axis  by  Bfs  =  tan'1(We/Ue).  and  then  rotating 
around  the  X  axis  by  -(90°-<p)  to  the  Xfs-Y(s-Zu  coordinate  system  for  each  radial  profile. 
In  this  coordinate  system,  for  each  data  point.  Y,5  is  the  distance  from  the  model  surface 
and  Zfs  =  0.0.  At  the  location  in  each  profile  where  the  resultant  velocity  Ur  =  Ure,  V  and 
W  are  zero  and  U  =  Ur=  Ure.  This  coordinate  axis  system  was  chosen  as  the  proper  local 
free-stream  coordinate  system  to  use  to  present  the  turbulence  measurements  and  derived 
quantities  so  that  they  have  meaning.  However,  because  the  measurements  were  made 
in  the  vertical  YT-7^r  plane,  any  distortions  in  the  locations  of  the  points  in  the  X- 
directions  were  ignored  after  the  data  were  transformed  to  the  local  free  stream  coordinate 
system.  Also,  when  finding  the  mean  velocity  gradients  in  the  Y(s-direction  for  the 
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derived  quantities  discussed  in  Section  4.3,  any  gradients  in  the  XT-direction  were 
ignored. 


As  a  basis  to  consider  the  validity  of  the  turbulent  measurement  results,  the 
realizability  conditions  outlined  in  Schumann  (1977)  were  utilized.  Basically,  the 
Reynolds  stress  tensor  for  each  data  point  must  satisfy  the  following: 

/j  *  u2  +  v2  +  w:  £  0 

/,  s  u2v2  -  (uv'.:  *  u2w2  -  (uvv)2 

-  \’2w2  -  (v»v)"  ^  0 
/,  ■  U2(v2U’2  -  Wi~) 

-  uv{uvw 2  -  VWUw) 

-  Uw(uVVW  -  v2uw)  ^  0 

If  any  of  the  above  conditions  are  not  satisfied  for  a  given  Reynolds  stress  tensor,  then 
the  validity  of  the  data  at  that  point  must  be  questioned.  The  present  measurements  w-ere 
tested  using  the  above  conditions,  and  for  all  the  data  at  all  the  points,  the  conditions 
were  satisfied.  Hence,  from  the  point  of  view  of  the  realizability  of  the  Reynolds  stress 
tensor,  the  data  are  valid. 

Figures  80  to  85  show  the  normal  stresses  nondimensionalized  by  Ure  for  all  the 
radial  profiles  at  X/L  =  0.7  and  X/L  =  0.8.  At  each  axial  location,  an  increase  in  these 
quantities  occurs  at  a  further  distance  from  the  surface,  given  by  Yfs.  as  <p  increases  from 
20°  to  70°.  At  ip  =  80°  and  90°,  this  distance  decreases.  This  corresponds  to  the  growth 
of  the  boundary  layer  circumferentially  towards  the  leeward  side,  as  seen  in  the  velocity 
profiles,  because  the  sudden  increase  in  each  quantity  occurs  at  the  edge  of  the  boundary 
layer.  The  turbulent  kinetic  energy,  given  by: 
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is  shown  in  Figures  86  and  87  for  both  axial  locations.  In  general,  the  turbulent  kinetic 
energy  increases  in  the  proximity  of  the  surface  in  the  same  manner  as  a  two-dimensional 
boundary'  layer,  and  there  appears  to  be  little  evidence  of  turbulent  energy  being 
convected  away  from  the  surface  by  the  vortical  flow. 

Figures  88  to  95  show  the  normal  stress  quantities  for  each  profile  for  X/L  =  0.7 
and  Figures  96  to  103  show  these  quantities  for  X/L  =  0.8.  For  both  axial  locations,  at 

<p  =  20°.  w2  is  greater  than  v:.  At  <p  =  30°.  these  quantities  are  approximately  equal  at 

X/L  =  0.7.  and  at  X/L  =  0.8,  v2  is  greater  than  w2  closer  to  the  surface.  At  ip=40°  and 

50°  for  both  axial  locations,  v2  is  greater  than  w2.  At  ip  =  60°.  v2  is  greater  thanw2 
until  close  to  the  surface,  where  these  quantities  are  equal.  At  ip  =  70°,  these  quantities 
are  approximately  equal  throughout  the  boundary  layer,  and  at  <p  =  80°  and  90°.  w2  is 

again  greater  than  v2 . 

The  locations  just  described  where  v2  is  greater  than  w2  occur  in  regions  in 
which  the  processes  of  separation  are  beginning  to  occur  along  the  surface  of  the  body 
and  vortical  flow  becomes  evident.  Although  Klebanoff  (1954)  showed  that  w2  is  greater 

than  v2  throughout  the  boundary  layer  for  turbulent  flow  over  a  flat  plate,  the  present 
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data  show  distortion  from  this  due  to  the  separation  phenomena  and  the  presence  of  three- 
dimensional  vortical  flow,  and  these  distortions  are  consistent  when  comparing  data  at 

both  axial  locations.  The  variations  of  the  relative  magnitudes  of  v2  and  w2  can  possibly 

be  attributed  to  the  folding  over  of  the  separating  flow.  Overall,  u 2  is  greater  than  eitherv2 

or  w2  in  all  the  profiles,  and  this  quantity  makes  up  most  of  the  turbulent  kinetic  energy 
through  the  boundary  layer. 

Shear  stress  quantities  -u\\  -uw,  and  -w  are  shown  in  Figures  104  to  109  for 

all  the  profiles  at  X/L=0.7  and  X/L=0.8.  As  with  the  normal  stresses,  the  shear  stress 
magnitudes  begin  to  increase  at  the  edge  of  the  boundary  layer. 

The  shear  stress  quantities  for  each  profile  are  shown  in  Figures  110  to  117  for 
X/L  =  0.7  and  in  Figures  118  to  125  for  X/L=0.8.  In  all  cases  except  X/L=0.8  and 

(p  =  90o,  uv  has  the  largest  magnitude,  but  the  behavior  of  uw  and  vw  are  still 
significant.  As  was  previously  described  concerning  the  relative  magnitude  of  w2  andv2 
at  the  same  <p  for  both  axial  locations,  the  behavior  of  uw  and  vw  is  consistent  when 
comparing  values  at  the  same  <p  at  both  axial  locations.  At  ip  =  20°,  vw  is  negative  at  the 
edge  of  the  boundary  layer  and  becomes  positive  at  approximately  Yfs/599=0.5,  anduw 
is  initially  negative  and  starts  to  become  positive.  At  <p  =  30°,  vw  is  positive  and  uw  has 
a  smaller  magnitude  than  vw.  At  tp  =  40°,  both  values  are  positive.  However,  at 
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<p  =  50°,  w  is  negative,  and  uw  and  w  have  approximately  the  same  magnitude.  At 
ip  =  60°  and  ip  =  70°,  uw  becomes  negative  below  YJ&99=0. 4.  At  <p  =  80°.  both 

quantities  are  negative  except  below  Yfs/6W=0.2,  where  vw  becomes  positive. 

Kreplin  and  Vollmers  (1981)  obtained  Reynolds  stress  components  on  a  6:1 
prolate  spheroid  at  X/L  =  0.64  and  a=  10°.  These  measurements  were  taken  in  regions 
of  converging  wall  streamlines  for  Re  =  7.2  x  106.  at  ip  =  60°  and  ip=45°.  Although  these 
data  were  obtained  at  a  higher  Reynolds  number,  comparisons  can  still  be  made  with  the 
present  data  for  regions  of  converging  wall  streamlines.  Based  on  the  discussion  in 
Section  4.1  concerning  results  of  oil  flows  done  by  Ahn  (1990),  these  regions  occur  at 
approximately  <p  =  30°  for  X/L  =  0.7  and  <p  =  20°  for  X/L=0.8. 

Figure  126  shows  the  comparison  of  the  u2  Reynolds  stresses.  The  data  at 

ip  =  20°  and  X/L=0.8  show  good  agreement  with  the  data  from  Kreplin  and  Vollmers. 
However,  the  data  at  <p  =  30°  and  X/L=0.7  give  lower  values.  As  seen  in  Figure  127. 

the  present  data  generally  give  higher  values  of  v2  than  the  data  from  Kreplin  and 

Vollmers.  The  comparison  of  w2  in  Figure  128  shows  that  the  present  data  at  X/L=0.7 

and  <p  =  30°  agree  well  with  the  values  from  Kreplin  and  Vollmers,  but  the  data  at 
X/L  =  0.8  and  tp  =  20°  show  higher  values. 

The  comparison  of  uv  in  Figure  129  shows  that  the  present  data  exhibit  the  same 
trends  as  the  data  from  Kreplin  and  Vollmers;  however,  the  magnitudes  of  the  present 


data  are  generally  higher.  The  comparison  of  uw  in  Figure  130  does  not  show  good 
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correlation  between  the  data.  The  present  data  are  generally  positive,  whereas  the  data 
from  Kreplin  and  Vollmers  are  negative.  The  comparison  of  w  in  Figure  131  shows 

somewhat  better  agreement,  although  the  values  from  Kreplin  and  Vollmers  have  a  larger 
magnitude  overall. 

Some  of  the  differences  in  the  previously  described  comparisons  can  be  attributed 
to  both  the  different  Reynolds  numbers  and  the  different  locations.  Also,  the  x-wire  data 
reduction  method  used  by  Kreplin  and  Vollmers  (1981)  neglected  tangential  cooling  and 
assumed  a  cosine  cooling  law  for  the  wires. 


4.3  Derived  Quantities 


The  Reynolds  stress  quantities  obtained  in  the  measurements  were  used  to  estimate 

various  derived  quantities,  defined  in  the  following 

discussion,  in  local  free  stream 

coordinates.  These  quantities  are  listed  below,  along 

with  the  estimated  uncertainties  of 

each  quantity: 

Turbulent  shear  stress  r 

±0.17Tmax 

Townsend's  structural  parameter  A1 

±0.1 8  A 1  max 

Kinematic  eddy  viscosity  vx 

Kinematic  eddy  viscosity  v, 

±0.37vz.m„ 

Anisotropic  constant  N 

±0.40Nrnax 

Mixing  length  Lm 

±0. 13Lm.max 
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Flow  gradient  angle  ag  ±0.095as.max 

Shear  stress  angle  a,  ±0.35asm<iX 

The  given  estimated  uncertainties  are  valid  for  data  within  the  inner  two-thirds  of  the 
boundary  layer.  At  the  edge  of  the  layer  and  outside  the  layer,  the  results  have 
increasing  uncertainties  due  to  the  small  magnitudes  of  the  velocity  gradients  and 
Reynolds  stress  quantities  involved  in  the  calculations.  Therefore,  values  of  the  derived 
quantities  at  Yfs/<5,^  greater  than  approximately  0.9  are  not  presented. 

The  total  turbulent  shear  stress  r  in  planes  parallel  to  the  Xh-Zti  plane  of  the  local 

free  stream  coordinate  system  can  be  obtained  from  uv  and  vw  by: 

=  y  (uvf  *  (vvvp 


The  results  for  X/L^O.7  and  X/L  =  0.8.  nondimensionalized  by  Ure:,  are  shown  in  Figure 
132  and  133.  The  magnitude  of  r  increases  through  the  boundary  layer  and  is  highest 
near  the  wall. 

Townsend's  structural  parameter  Al  relates  the  turbulent  kinetic  energy  k  to  the 
Reynolds  shear  stresses,  and  is  given  by  the  ratio  of  t  to  k  as  follows: 


Al 


\l{uv)2  +  (w|2 
(u*  +  V2  *  w2) 


This  parameter  has  been  found  to  be  nearly  constant  at  a  value  of  0.15  in  two- 
dimensional  boundary  layers. 
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The  results  for  this  three-dimensional  boundary  layer  at  X/L=0.7  and  X/L  =  0.8 
are  shown  in  Figures  134  and  135.  At  the  edges  of  the  layer.  A1  tends  to  vary  from 
approximately  0.05  to  0.17.  However,  as  Yfs/<5^  decreases,  the  values  for  A1  remain 
within  a  band  from  0.13  to  0.17.  Thus,  although  the  parameter  is  distorted  somewhat 
from  the  value  of  0.15  for  two-dimensional  layers,  the  ability  of  turbulence  to  generate 
transport  of  momentum  in  this  three-dimensional  layer  appears  to  be  similar  to  that  of  a 
two-dimensional  layer,  and  assuming  a  value  of  0. 15  for  Townsend's  structural  parameter 
A1  is  not  wholly  inaccurate. 

The  kinematic  eddy  viscosities  in  the  x  direction  and  in  the  z  direction  for  the 
local  free  stream  coordinate  system  are  given  by: 

-uv 

V  =  - 

'dU' 

k  &y  / 

-vw 

V.  =  - 

■  IdW' 

1  dy, 


The  gradients  in  the  Yfs  direction  were  estimated  by  fitting  a  parabola  through  five 
consecutive  points  and  finding  the  derivative  at  the  third  (middle)  point. 

The  results  for  vx.  shown  in  Figure  136  and  137.  are  everywhere  positive  in  the 
boundary  layer.  Above  about  Yf,/<SW  =  0.3.  the  x-eddy  viscosity  appears  to  be  varying 
non-linearly  with  respect  to  Iog(Yfs/6w).  Below  this,  the  variation  is  more  linear.  The 
results  for  v,,  shown  in  Figures  138  and  139.  are  more  difficult  to  interpret.  In  some 
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cases  a  linear  variation  with  respect  to  log(Yfs/<5w)  is  evident;  however,  there  is  no 
overall  consistency  in  the  behavior.  Most  of  the  values  are  in  a  band  within  ±0.004 
nr/s. 

To  determine  whether  a  correlation  exists  between  the  eddy  viscosity  and  the 
displacement  thickness,  the  x  eddy  viscosity  was  normalized  as  vx/(Ure6‘).  According 
to  the  Clauser  eddy  viscosity  model  for  a  flat-plate  zero  pressure  gradient  flow,  as  taken 
from  Schetz  (1984).  this  normalization  leads  to  a  value  of  0.018.  However,  for  the 
present  flow,  this  normalization  does  not  appear  to  provide  any  collapse  of  the  curves, 
as  shown  in  Figures  140  to  141 .  Within  the  boundary  layer,  the  values  for  X/L=0.7  are 
generally  between  0.010  and  0.020,  and  the  values  for  X/L=0.8  are  between  0.005  and 
0.010.  Generally,  as  <p  increases  from  20°  to  50°,  the  values  decrease.  As  <p  increases 
further  to  90°,  the  values  increase. 

Isotropic  eddy  viscosity  models  have  been  used  by  researchers  to  relate  the 
transfer  of  momentum  by  velocity  fluctuations  to  local  mean  flow  gradients.  Such  models 
utilize  an  eddy  viscosity  in  a  formulation  analogous  to  the  kinematic  viscosity  in  laminar 
flow.  However,  the  eddy  viscosity  is  not  a  thermophysical  property  of  the  fluid,  but  is 
a  property  of  the  flow.  In  order  to  test  the  validity  of  the  isotropic  assumption  as  applied 
to  the  present  three-dimensional  turbulent  flow,  the  anisotropic  constant  N  =  vz/vx  can  be 
examined.  If  the  isotropic  assumption  is  valid,  then  N=l.  How-ever,  as  shown  in 
Figures  142  and  143,  the  anisotropic  constant  varies  through  the  boundary  layer  in  a  band 
within  ±2.  Thus,  the  eddy  viscosity  appears  not  to  be  isotropic  in  this  turbulent  three- 
dimensional  boundary  layer. 
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In  addition  to  eddy  viscosity  methods,  the  mixing  length  Lm  is  used  to  relate 
turbulent  shear  to  mean  velocity  gradients.  The  mixing  length  is  defined  by  some 
interaction  distance  between  eddies  analogous  to  a  molecular  mean  free  path,  and  for  a 
three-dimensional  boundary  layer  L,,,  is  given  by: 


For  the  outer  region  of  the  boundary  layer,  where  most  of  the  present  data  were  obtained, 
a  generally-accepted  expression  for  the  mixing  length  for  two-dimensional  zero  pressure 
gradient  boundary'  layers  is  given  by  Lm;D  =  0.096  (Schetz.  1984),  resulting  in  a  constant 
mixing  length  for  such  a  boundary  layer.  To  determine  the  applicability  of  this 
expression  to  the  present  data,  the  ratio  of  the  mixing  length  L,,,  to  the  two-dimensional 
mixing  length  Lra:D  was  found  and  the  results  are  shown  in  Figure  144  and  145.  In 
general,  the  expression  for  two-dimensional  boundary  layers  overestimates  the  value  of 
the  mixing  length,  and  the  ratio  ranges  from  0.5  to  1.0.  However,  the  assumption  of  a 
constant  mixing  length  for  the  outer  region  of  the  layer  appears  to  be  valid. 

The  flow  gradient  angle  ag  and  the  shear  stress  angle  a5,  given  by: 

'  dW' 

a  =  arctan-— — 

*  (— ) 

a  =  arctan 
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were  calculated  in  an  attempt  to  determine  if  the  mean  velocity  gradients  align  in  the 
same  directions  as  the  shear  stresses.  The  results  for  ag  and  as  are  given  in  Figures  146 
to  149  for  X/L  =  0.7  and  X/L  =  0.8.  The  flow  gradient  angle  tends  to  increase  through 
the  layer,  except  at  tp  =  60°,  70°.  and  80°,  whereas  the  shear  angle  tends  to  decrease. 
Between  <p  =  50°  and  <p  =  80°,  the  trends  in  both  ag  and  as  appear  consistent  for  a  given 
profile;  however,  the  actual  magnitudes  differ  greatly.  Overall,  a  comparison  between 
the  flow  gradient  angle  and  the  shear  stress  angle  does  not  yield  consistent  results,  except 
that  the  directions  of  the  mean  velocity  gradients  and  the  corresponding  shear  stresses  do 
not  appear  aligned. 

In  an  effort  to  examine  the  lag  or  lead  of  the  shear  stress  angle  with  respect  to  the 
mean  velocity  flow  gradient  angle,  the  values  of  the  difference  ag-as  were  found  and  are 
shown  in  Figures  150  and  151.  Overall,  the  values  generally  vary  from  negative  at  the 
edge  of  the  boundary  layer  to  positive  at  the  surface,  indicating  that  at  the  edge  of  the 
layer  the  shear  stress  angle  leads  with  respect  to  the  flow  gradient  angle  while  at  the 
surface  the  shear  stress  angle  lags  when  compared  to  the  flow  gradient  angle.  A  notable 
exception  appears  at  ip  =  80°  for  both  axial  locations,  in  which  cases  the  values  tend  to 
decrease  through  the  layer  and  then  increase  close  to  the  surface.  At  each  axial  location, 
the  curves  for  <p  =  50°  and  <p  =  60°  have  the  least  variation  through  the  boundary  layer. 
As  shown  in  Figure  75,  these  are  regions  where  the  displacement  thickness  is  largest. 
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The  mean  velocity  field  and  turbulence  measurements  were  made  along  the 
leeward  side  of  a  6:1  prolate  spheroid,  in  regions  leading  to  and  including  separation. 

For  the  cases  of  Re=1.3xl06  at  a=10°  and  15°,  separation  results  in  weak 
vortical  flow,  as  shown  by  the  magnitudes  of  the  secondary  velocities  and  by  the 
maximum  vorticity  nx/U„  of  0.10/cm  in  the  center  of  the  region  of  vortical  flow  at 
X/L  =  0.9.  Based  upon  results  of  surface  skin  friction  measurements  from  Kreplin,  et  al. 
(1980).  the  flow  in  each  case  may  possibly  be  laminar  up  to  separation. 

For  Re  =  4.0 x  106,  transition  occurs  along  the  forebody,  and  the  flow  is  turbulent 
in  the  locations  of  the  measurements.  At  cr  =  10°,  separation  begins  to  occur  between 
X/L  =  0.8  and  X/L=0.9  and  the  resulting  vortical  flow  is  stronger  than  in  the  lower 
Reynolds  number  case,  with  larger  secondary  velocities  and  a  maximum  vorticity  nx/U„ 
greater  than  0.16/cm  at  X/L  =  0.9.  At  a=  15°,  separation  occurs  further  forward  and  the 
resultant  vortical  flow  is  the  strongest  of  any  of  the  cases,  with  a  maximum  vorticity 
nx/U„  at  X/L  =  0.9  of  0.30/cm.  At  both  angles  of  attack,  the  diffusion  of  vorticity  away 
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from  the  body  is  contained  by  the  presence  of  the  vortical  flow.  The  centers  of  the 
regions  of  strong  vortical  flow  appear  to  move  away  from  the  surface  slightly  by 
X/L  =  0.9. 


Estimates  of  the  boundary  layer  thickness  <5^  and  the  displacement  thickness  <5* 
indicate  the  growth  of  the  boundary  layer  circumferentially  and  axially  along  the  body 
for  q=  10°  and  Re  =  4.0x  106.  Comparisons  with  data  from  Meier,  et  al.  (1984a)  show 
that  at  a  higher  Reynolds  number,  the  displacement  thickness  increases  and  the  location 
of  maximum  displacement  thickness  shifts  toward  the  leeward  side.  The  locations  of 
negative  crossflow  W/Ure  (in  the  local  body  coordinate  system)  indicated  by  a  change  in 
the  sign  of  W/Ure  occur  further  away  from  the  leeward  side  at  a  higher  Reynolds  number. 
This  flow  reversal  is  coincident  with  a  large  increase  in  the  displacement  thickness. 

Skewing  of  the  boundary  layer  is  indicated  by  the  variation  of  the  direction  of  the 
velocity  vectors  through  the  layer.  Maximum  skewing  occurs  at  <p  =  60",  independent  of 
axial  location  and  Reynolds  number,  and  this  location  is  closer  to  the  leeward  side  than 
the  location  of  maximum  displacement  thickness.  Again,  comparison  with  Meier,  et  al. 
(1984a)  shows  good  agreement,  considering  Reynolds  number  effects. 

Turbulence  measurements  in  the  regions  leading  up  to  separation  show  that  the 
distribution  of  Reynolds  stresses  is  distorted  compared  to  the  distribution  of  a  two- 
dimensional  boundary  layer.  This  distortion  can  be  attributed  to  the  separation 
phenomena  and  the  presence  of  three-dimensional  vortical  flow.  Comparing  the 
distributions  at  X/L  =  0.7  with  those  at  X/L  =  0.8  indicates  that  they  are  consistent  axially 
along  the  body.  Comparison  of  Reynolds  stresses  in  regions  of  converging  streamlines 
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with  those  from  Kreplin  and  Vollmers  (1981)  show  good  agreement  in  most  cases  for  the 
Reynolds  normal  stresses;  however,  the  correlation  is  not  as  good  for  the  Reynolds  shear 
stresses. 

Estimates  of  Townsend's  structural  parameter  A 1  relating  turbulent  kinetic  energy 
to  the  Reynolds  shear  stresses  indicates  that  the  distribution  is  similar  to  that  of  a  two- 
dimensional  layer  with  A1  =0.15.  However,  estimates  of  eddy  viscosity  in  the  x  and  z 
directions  and  the  resulting  anisotropy  constant  N  indicate  that  eddy  viscosity  is  not 
isotropic  through  the  boundary  layer.  Hence,  isotropic  models  for  eddy  viscosity  will  not 
give  accurate  results.  Estimates  of  the  mixing  length  compared  to  a  generally-accepted 
theoretical  two-dimensional  model  show  that  the  model  overestimates  the  mixing  length 
for  the  outer  region.  A  comparison  of  the  flow  gradient  angles  and  the  shear  stress 
angles  shows  that  the  mean  velocity  gradients  and  the  shear  stresses  are  not  aligned  in  this 
three-dimensional  boundary  layer.  Overall,  the  shear  stress  angle  leads  the  mean  velocity 
flow  gradient  angle  at  the  edge  of  the  boundary  layer,  while  at  the  surface  the  shear  stress 
angle  lags  compared  to  the  flow  gradient  angle. 

For  this  study,  turbulence  measurements  were  obtained  in  the  outer  region  of  the 
boundary  layer  only,  using  an  x-wire  mounted  on  a  traverse  external  to  the  model. 
Future  work  could  include  measurements  closer  to  the  surface,  using  hot-wire 
anemometry  and  laser  doppler  velocimetry  with  probes  mounted  and  traversed  inside  the 
model.  The  present  data  could  be  used  for  comparison  purposes.  Since  this  is  part  of 
an  overall  study  of  transition  and  separation,  future  measurements  should  be  obtained  in 
both  steady  and  unsteady  flow. 
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In  addition  to  turbulence  measurements,  skin  friction  measurements  are  important 
in  the  study  of  transition  and  separation.  Such  data  would  confirm  results  from  surface 
oil  flow  measurements  and  would  aid  in  the  analysis  of  turbulence  data  for  the  inner 
region  and  lamina  -  sublayer  of  the  boundary  layer  in  both  steady  and  unsteady  cases. 
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Illustrations 


Illustrations 


Wind  Tunnel. 


Figure  2. 


Figure  3.  6:1  Prolate  Spheroid  Model  Mounted  in  Wind  Tunnel. 
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Figure  7.  X-Wire  Probe  Support  with  Pitch  and  Roil  Capabilities. 
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Figure  8.  Diagram  of  Instrumentation  for  X-Wire  Probe. 
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Figure  12. 


Results  of  Surface  Oil-Flo' 
Re=4.0x  106,  (Ahn.  1990). 


1 1  lustrations 


Visualization  Experiments.  a=10°. 
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Figure  13. 


Results  of  Surface  Oil-Flow  Visualization  Experiments. 
Re=4.0x  106.  (Ahn,  1990). 
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Figure  14. 


Results  of  Surface  Oil-Flow  Visualization  Experiments. 
Re=4.0x  106.  (Ahn.  1990). 
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Figure  15.  Diagram  of  Cross-Section  of  Separation  Sheet  of  Fluid 
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Figure  !  6.  Secondary  Velocity  Vectors  from  Yawhead  Probe  Measurements,  in  Tunnel 
Coordinates.  X/L  =  0.9,  a=15°.  Re=4.0xl0ft,  Grid  Spacmg=  1 ,27cm 
(0.5in).  Compare  with  Figure  32  to  see  Results  of  Rotation. 
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Figure  17. 
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Secondary  Velocity  Vectors  from  Yawhead  Probe  Measurements. 
X/L  =  0.6.  a=  10°,  Re=  1 .3  x  lQft.  Grid  Spacing  of  1 ,27cm  (0.5in). 
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Secondary  Velocity  Vectors  from  3awhead  Probe  Measurements. 
X/L=0.7,  a  =  10°,  Re=  1 .3  x  10ft.  Grid  Spacing  of  1 ,27cm  (0.5in). 
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f  igure  19.  Secondary  Velocity  Vectors  from  Yauhead  Probe  Measurements 
-^',L=0.8.  a=  10°.  Re=  1 .3 x  106.  Grid  Spacing  of  1 .27cm  (0.5in). 
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Figure  20.  Secondary  Velocity  Vectors  from  Yawhcad  Probe  Measurements. 
X/L  =  0.9,  cl  —  1 0 c .  Re=1.3xl06.  Grid  Spacing  of  1.27cm  (0.5in). 
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Secondary  Velocity  Vectors  from  Yaw  head  Probe  Measurements 
X/L  =  0.6.  a=15°.  Re=1.3xl06,  Grid  Spacing  of  1.27cm  (0.5in). 
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Figure  23 


B 


Secondary  Velocity  Vectors  from  Yawhead  Probe  Measurements. 
X/L  =  0.8,  a=  15°.  Re=  1 .3 x  10*.  Grid  Spacing  of  1 .27cm  tu.5in). 
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Figure  24.  Secondary  Velocity  Vectors  from  Yawhead  Probe  Measurements  Ttw 
X/L  =  0.9.  a=  15°.  Re=  1 .3 x  106.  Grid  Spacing  of  1 .27cm  (0.5in). 
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B 


- > 

0.33U. 


Secondary  Velocity  Vectors  from  Yawhead  Probe  Measurements, 
X/L  =  0.8.  a=  10°,  Re  =  4.0x  106.  Grid  Spacing  of  1 ,27cm  (0.5in). 
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Figure  28. 


Secondary  Velocity  Vectors  from  Yawhead  Probe  Measurements. 
X/L  =  0.9,  a=  10°.  Re=4.0x  106.  Grid  Spacing  of  1 .27cm  (0.5in). 
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Figure  29.  Secondary  Velocity  Vectors  from  Yawhead  Probe  Measurements. 
X/L  =  0.6,  a  =  15°,  Re=4.0x  106.  Grid  Spacing  of  1.27cm  (0.5in). 
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Secondary  Velocity  Vectors  from  Yawhead  Probe  Measurements, 
X/L  =  0.7.  a=  15°,  Re  =  4.0x  106,  Grid  Spacing  of  1 .27cm  (0.5in). 
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Figure  31 . 
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Secondary  Velocity  Vectors  from  Yawhead  Probe  Measurements. 
X/L  =  0.8.  a-  15°.  Re=4.0x  lO6.  Grid  Spacing  of  1 ,27cm  (0.5in). 
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Figure  32. 


Secondary  Velocity  Vectors  from  Yawhead  Probe  Measurements. 
X/L  =  0.9.  a=  15°,  Re  =  4.0x  106.  Grid  Spacing  of  1.27cm  (0.5in). 
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Figure  49. 


Locations  of  X-Wire  Measurement  Points  in  Horizontal  and  Vertical 


Profiles  and  Locations  of  Interpolated  Points  in  Radial  Profiles. 
X/L=0.7. 
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Fieure  50. 


Locations  of  X-Wire  Measurement  Points  in  Horizontal  and  Vertical 
Profiles  and  Locations  of  Interpolated  Points  in  Radial  Profiles. 
X/L=0.8. 
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Figure  51.  Comparison  of  Mean  Velocity  Components  in  Tunnel  Coordinates  from 
X-Wire  Measurements  and  Yawhead  Probe  Measurements,  ^  =  -3. 81cm. 
X/L  =  0.7,  Re=4.0xl06,  a  =  10°. 
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Figure  52. 


Comparison  of  Mean  Velocity  Components  in  Tunnel  Coordinates  from 
X-Wire  Measurements  and  Yawhead  Probe  Measurements,  7^  =  -5. 08cm, 
X/L=0.7,  Re=4.0xl06,  a=10°. 
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Figure  53.  Comparison  of  Mean  Velocity  Components  in  Tunnel  Coordinates  from 
X-Wire  Measurements  and  Yawhead  Probe  Measurements.  Zt  =  -10.  15cm. 
X/L  =  0.7,  Re  =  4.0xl06,  a  =  10°. 
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Figure  54.  Comparison  of  Mean  Velocity  Components  in  Tunnel  Coordinates  from 
X-Wire  Measurements  and  Yawhead  Probe  Measurements,  ZT  =  -1.27cm, 
X/L=0.8,  Re=4.0xl06,  ct  =  10o. 
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Comparison  of  Mean  Velocity  Components  in  Tunnel  Coordinates  from 
X-Wire  Measurements  and  Yawhead  Probe  Measurements,  ^=-3. 81cm, 
X/L  =  0.8,  Re  =  4.0xl06,  a  =  10°. 
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Figure  56.  Comparison  of  Mean  Velocity  Components  in  Tunnel  Coordinates  from 
X-Wire  Measurements  and  Yawhead  Probe  Measurements,  ^  =  -5. 08cm. 
X/L=0.8,  Re=4.0xl06,  a  =  10°. 


Illustrations 


108 


» 


» 


t 


Velocity /  Ure 


Figure  57.  Comparison  of  Mean  Velocity  Components  in  Tunnel  Coordinates  from 
X-Wire  Measurements  and  Yawhead  Probe  Measurements,  Zr  =  -6.35cm, 
X/L=0.8,  Re=4.0x  106,  a=10°. 
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Figure  58.  Local  Body  (Xb-Yb-Zb)  Coordinate  System,  with  Xb  Parallel 
Longitudinal  Axis  of  Model  and  Yb  out  from  Surface. 
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Figure  59. 
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X-Wire  Velocity  Profiles  and  U-W  Crossflow  Vectors,  <p  =  20°,  X/L=0.7, 
Re=4.0xl06,  a  =  lfl°.  Local  Body  Coordinates,  Direction  of  Arrow  on 
Hodograph  Plot  Indicates  Variation  from  Wall  to  Free  Stream. 
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Figure  60.  X-Wire  Velocity  Profiles  and  U-W  Crossflow  Vectors,  ip  =  30°,  X/L=0.7, 
Re =4. Oxl 06,  a=10°,  Local  Body  Coordinates,  Direction  of  Arrow  on 
Hodograph  Plot  Indicates  Variation  from  Wall  to  Free  Stream. 
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Figure  61 . 
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X-Wire  Velocity  Profiles  and  U-W  Crossflow  Vectors,  ip=40°,  X/L=0.7, 
Re=4.0xl06,  a  =  10°,  Local  Body  Coordinates,  Direction  of  Arrow  on 
Hodograph  Plot  Indicates  Variation  from  Wall  to  Free  Stream. 
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Figure  62. 
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X-Wire  Velocity  Profiles  and  U-W  Crossflow  Vectors.  <p  =  50°,  X/L=0.7, 
Re  =  4.0xl06,  a=10°.  Local  Body  Coordinates,  Direction  of  Arrow  on 
Hodograph  Plot  Indicates  Variation  from  Wall  to  Free  Stream. 
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Figure  63.  X-Wire  Velocity  Profiles  and  U-W  Crossflow  Vectors,  <p  =  6  \  X/L  =  0.7, 
Re  =  4.0xl06,  a  =  10°,  Local  Body  Coordinates,  Direction  of  Arrow  on 
Hodograph  Plot  Indicates  Variation  from  Wall  to  Free  Su.^m. 
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Figure  64.  X-Wire  Velocity  Profiles  and  U-W  Crossflow  Vectors,  tp=70°,  X/L=0.7, 
Re=4.0xl06,  a  =  10°,  Local  Body  Coordinates,  Direction  of  Arrow  on 
Hodograph  Plot  Indicates  Variation  from  Wall  to  Free  stream. 
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Figure  65.  X-Wire  Velocity  Profiles  and  U-W  Crossflow  Vectors,  <p  =  80°,  X/L=0.7, 
Re=4.0xl06,  a=10°,  Local  Body  Coordinates,  Direction  of  Arrow  on 
Hodograph  Plot  Indicates  Variation  from  Wall  to  Free  Stream. 
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Figure  66.  X-Wire  Velocity  Profiles  and  U-W  Crossflow  Vectors,  <p  =  90°,  X/L=0.7, 
Re=4.0xl06.  a  =  10°,  Local  Body  Coordinates,  Direction  of  Arrow  on 
Hodograph  Plot  Indicates  Variation  from  Wall  to  Free  Stream. 
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Figure  67.  X-Wire  Velocity  Profiles  and  U-W  Crossflow  Vectors,  ip  =  20°,  X/L=0.8, 
Re  =  4.0xl06,  a=10°,  Local  Body  Coordinates,  Direction  of  Arrow  on 
Hodograph  Plot  Indicates  Variation  from  Wall  to  Free  Stream. 
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Figure  68.  X-Wire  Velocity  Profiles  and  U-W  Crossflow  Vectors,  <p  =  30°,  X/L=0.8, 

Re=4.0xl06,  a=10°,  Local  Body  Coordinates,  Direction  of  Arrow  on 
Hodograph  Plot  Indicates  Variation  from  Wall  to  Free  Stream. 
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Figure  69.  X-Wire  Velocity  Profiles  and  U-W  Crossflow  Vectors,  <p  =  40°,  X/L=0.8, 
Re  =  4.0xl06,  a  =  10°,  Local  Body  Coordinates,  Direction  of  Arrow  on 
Hodograph  Plot  Indicates  Variation  from  Wall  to  Free  Stream. 
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Figure  70.  X-Wire  Velocity  Profiles  and  U-W  Crossflow  Vectors,  (p  =  50°,  X/L=0.8, 
Re  =  4.0xl06,  a  =  10°,  Local  Body  Coordinates,  Direction  of  Arrow  on 
Hodograph  Plot  Indicates  Variation  from  Wall  to  Free  Stream. 
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Figure  71.  X-Wire  Velocity  Profiles  and  U-W  Crossflow  Vectors,  tp  =  60°,  X/L=0.8, 
Re=4.0xl06,  a  =  10°,  Local  Body  Coordinates,  Direction  of  Arrow  on 
Hodograph  Plot  Indicates  Variation  from  Wall  to  Free  Stream. 
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Figure  72.  X-Wire  Velocity  Profiles  and  U-W  Crossflow  Vectors,  <p  =  70°,  X/L=0.8, 
Re  =  4.0xl06,  a=10°,  Lx>cal  Body  Coordinates,  Direction  of  Arrow  on 
Hodograph  Plot  Indicates  Variation  from  Wall  to  Free  Stream. 
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Figure  73. 
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X-Wire  Velocity  Profiles  and  U-W  Crossflow  Vectors,  <p  =  80°,  X/L=0.8, 
Rr  =4.0xl06,  a  =  10°,  Local  Body  Coordinates,  Direction  of  Arrow  on 
Hodograph  Plot  Indicates  Variation  from  Wall  to  Free  Stream. 
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Figure  74.  X-Wire  Velocity  Profiles  and  U-W  Crossflow  Vectors.  <p  =  90°,  X/L=0.8. 

Re=4.0xl06,  a=  10°,  Local  Body  Coordinates,  Direction  of  Arrow  on 
Hodograph  Plot  Indicates  Variation  from  Wall  to  Free  Stream. 
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Figure  75.  Displacement  Thickness  S'  with  respect  to  Circumferential  Angle  <p. 

Re  =  4.0xl06.  a=10°,  Compared  to  DFVLR  Data  at  Re  =  7.2xl06 
(Meier,  ei  al.  ,1984a). 


Illustrations 


127 


Yb  (cm) 


-10.0  -5.0  0.0  5.0  10.0  15.0  20.0 

P  (degrees) 


Figure  76.  Profiles  of  Flow  Angle  B  =  arctan(W/U),  X/L  =  0.7,  Re =4.  Ox  106,  a=  10°, 

Local  Body  Coordinates. 
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Figure  77.  Profilesof  Flow  Angle  J3  =  arctan(W/U),  X/L=0.8.  Re=4.0x  106.  a=  10°, 
Local  Body  Coordinates. 


Illustrations 


129 


A  (3  ( degrees ) 


0.0  20.0  40.0  60.0  80.0  100.0 

({>  ( degrees ) 


Figure  78.  Difference  Between  Flow  Angle  of  Local  Free-Stream  Velocity  and  Flow 
Angle  of  Velocity  at  the  Surface  with  respect  to  Circumferential  Angle  »p, 
Re=4.0xl06,  a  =  10°,  Compared  to  DFVLR  Data  at  Re  =  7.2xl06 
(Meier,  et  al. ,  1984a). 
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Figure  79.  Local  Free-Stream  (Xfs-Yfj-Zfs)  Coordinate  System,  with  Xfs  Aligned  with 
the  Local  Free-Stream  Velocity  and  Yfs  out  from  Surface  of  Model. 
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Figure  80. 
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Profiles  of  u2  Reynolds  Stress,  X/L=0. 7,  Re  =  4. Ox  106,  a  =  10°,  Local 
Free- Stream  Coordinates. 
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Figure  81. 


Profiles  of  v2  Reynolds  Stress.  X/L=0. 7,  Re=4. Ox  1 06,  a 
Free-Stream  Coordinates. 
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Figure  83.  Profilesof  u 2  Reynolds  Stress,  X/L  =  0. 8,  Re=4. Ox  106,  a=  10°,  Local 
Free -Stream  Coordinates. 
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Figure  84.  Profilesof  v2  Reynolds  Stress,  X/L  =  0. 8,  Re=4. Ox  106,  a  =  10°,  Local 
Free-Stream  Coordinates. 
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Figure  85.  Profilesof  w 2  Reynolds  Stress,  X/L=0. 8,  Re=4. Ox  106,  a=  10°,  Local 
Free-Stream  Coordinates. 
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Profiles  of  Turbulent  Kinetic  Energy  k  with  respect  to  X/L=0.8, 

Re=4.0xl06,  a=10°,  Local  Free-Stream  Coordinates. 
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Figure  88.  Profiles  of  Reynolds  Normal  Stresses  with  respect  to  YJS,#,  <p  = 
X/L=0.7,  Re=4.0xl06,  a  =  10°,  Local  Free-Stream  Coordinates. 
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Figure  89.  Profiles  of  Reynolds  Normal  Stresses  with  respect  to  y- 

X/L=0.7,  Re  =  4.0xl06,  a=10°,  Local  Free-Stream  Coordinates. 
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Figure  90.  Profiles  of  Reynolds  Normal  Stresses  with  respect  to  Yfs/ <5^,  <p=40 
X/L=0.7,  Re=4.0xl06,  a  =  10°,  Local  Free-Stream  Coordinates. 
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Figure  91.  Profiles  of  Reynolds  Normal  Stresses  with  respect  to  <p  =  50°, 

X/L=0.7,  Re=4.0xl06,  a=10°,  Local  Free-Stream  Coordinates. 
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Figure  92.  Profiles  of  Reynolds  Norma!  Stresses  with  respect  to  Y <p =  60° 
X/L  =  0.7,  Re  =  4.0xl06,  a=10°,  Local  Free-Stream  Coordinates. 


Illustrations 


Figure  93. 


Profiles  of  Reynolds  Normal  Stresses  with  respect  to  <p  = 

X/L=0.7,  Re  =  4.0xl06,  a=10°.  Local  Free-Stream  Coordinates. 
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Figure  94.  Profiles  of  Reynolds  Normal  Stresses  with  respect  to  YJS^,  <p 
X/L  =  0.7,  Re=4.0xl06,  a=10°.  Local  Free-Stream  Coordinates 
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Figure  95.  Profiles  of  Reynolds  Normal  Stresses  with  respect  to  Yh/6<#,  t p- 
X/L=0.7,  Re  =  4.0xl06,  a=10°,  Local  Free-Stream  Coordinates. 
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Figure  96.  Profiles  of  Reynolds  Normal  Stresses  with  respect  to  <p 

X/L=0.8,  Re=4.0xl06,  a  =  10°.  Local  Free-Stream  Coordinates 
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Figure  97.  Profiles  of  Reynolds  Normal  Stresses  with  respect  to  YJS^,  ip  =  3 
X/L=0.8,  Re=4.0xl06,  a=10°.  Local  Free-Stream  Coordinates. 
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Figure  98.  Profiles  of  Reynolds  Normal  Stresses  with  respect  to  <p: 


X/L=0.8,  Re  =  4.0xl06,  a=10°.  Local  Free-Stream  Coordinates. 
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Figure  99.  Profiles  of  Reynolds  Normal  Stresses  with  respect  to  <p: 

X/L=0.8,  Re=4.0xl06,  a  =  10°,  Local  Free -Stream  Coordinates. 
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Figure  100.  Profiles 
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of  Reynolds  Normal  Stresses  with  respect  to  Y J&w  V 
8,  Re=4.0xl06,  a  =  10°,  Local  Free-Stream  Coordinates. 
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Figure  101.  Profiles  of  Reynolds  Normal  Stresses  with  respect  to  YfJS^, 

X/L  =  0.8,  Re=4.0xl06,  a=10°,  Local  Free-Stream  Coordinates 
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Figure  103.  Profiles  of  Reynolds  Normal  Stresses  with  respect  to  YJS^,  <p 
X/L  =  0.8,  Re=4.0xl06,  a=10°,  Local  Free-Stream  Coordinates 
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Figure  104.  Profiles  of  uv  Reynolds  Stresses,  X/L=0.7,  Re  =  4.0xl06,  a  = 
Local  Free-Stream  Coordinates. 
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Figure  105.  Profiles  of  uw  Reynolds  Stresses,  X/L  =  0.7,  Re=4.0xl06,  a  =  10° 
Local  Free-Stream  Coordinates. 
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Figure  106  Profiles  of  wv  Reynolds  Stresses,  X/L  =  0.7,  Re  =  4.0xl0h,  a=l 
Local  Free-Stream  Coordinates. 
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Figure  107.  Profiles  of  uv  Reynolds  Stresses.  X/L  =  0.8.  Re  =  4.0xl06,  a=10 
Local  Free-Stream  Coordinates. 
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Figure  108.  Profiles  of  uw  Reynolds  Stresses.  X/L  =  0.8.  Re  =  4.0xl06.  a=10°. 
Local  Free-Stream  Coordinates. 
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Figure  109.  Profiles  of  vw  Reynolds  Stresses.  X/L  =  0.8.  Re=4.0xl06,  a=10°, 
Local  Free-Stream  Coordinates. 
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Figure  110.  Profiles  of  Reynolds  Shear  Stresses  with  respect  to  Yf <P 
X/L  =  0.7,  Re  =  4.0xl06,  a  =  10°,  Local  Free-Stream  Coordinates 
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Figure  111.  Profiles  of  Reynolds  Shear  Stresses  with  respect  to  Y < p 
X/L=0.7,  Re=4.0xl0ft,  q  =  10°.  Local  Free-Stream  Coordinates 
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Profiles  of  Reynolds  Shear  Stresses  with  respect  to  <p=40°, 

X/L=0.7,  Re  =  4.0xl06,  a=10°,  Local  Free-Stream  Coordinates. 
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Figure  113.  Profiles  of  Reynolds  Shear  Stresses  with  respect  to  Y <p 
X/L=0.7,  Re=4.0xl06,  a  =  10°,  Local  Free-Stream  Coordinates 
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Figure  114.  Profiles  of  Reynolds  Shear  Stresses  with  respect  to  YJS^,  <p 
X/L=0.7,  Re=4.0xl06,  a  =  10°,  Local  Free-Stream  Coordinates 
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Figure  115.  Profiles  of  Reynolds  Shear  Stresses  with  respect  to  Y fs/<5^,  <p: 

X/L=0.7,  Re  =  4.0xl06,  cr  =  10°,  Local  Free-Stream  Coordinates. 
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Figure  116.  Profiles  of  Reynolds  Shear  Stresses  with  respect  to  <p  =  80°, 

X/L  =  0.7,  Re  =  4.0xl06,  a=10°,  Local  Free-Stream  Coordinates. 
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Figure  117.  Profiles  of  Reynolds  Shear  Stresses  with  respect  to 

X/L  =  0.7,  Re=4.0xl06,  a=10°,  Local  Free-Stream  Coordinates 
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Figure  118.  Profiles  of  Reynolds  Shear  Stresses  with  respect  to  Yf JS^. 

X/L  =  0.8,  Re  =  4.0xl06,  a=10°,  Local  Free-Stream  Coordinates 
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Figure  119.  Profiles  of  Reynolds  Shear  Stresses  with  respect  to  Yfs/<5W,  <p 
X/L  =  0.8,  Re  =  4.0xl06,  a=10°,  Local  Free-Stream  Coordinates 
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Figure  120.  Profiles  of  Reynolds  Shear  Stresses  with  respect  to  YJS^,  <p=40 
X/L=0.8,  Re=4.0xl06,  a=10°,  Local  Free-Stream  Coordinates. 
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Figure  121.  Profiles  of  Reynolds  Shear  Stresses  with  respect  to  YJS^,  tp 
X/L  =  0.8,  Re  =  4.0xl06,  a=10°,  Local  Free-Stream  Coordinates. 
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Figure  122.  Profiles  of  Reynolds  Shear  Stresses  with  respect  to  Y JSgg,  ip  =  60‘ 
X/L=0.8,  Re=4.0xl06,  a  =  10°,  Local  Free-Stream  Coordinates. 
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Figure  123.  Profiles  of  Reynolds  Shear  Stresses  with  respect  to  <*> 

X/L=0.8,  Re=4.0xl06,  a=10°,  Local  Free-Stream  Coordinates 
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Figure  124.  Profiles  of  Reynolds  Shear  Stresses  with  respect  to  99'  *P: 

X/L=0.S,  Rc=4.0xl0*,  a  =  10°,  Local  Free  Stream  Coordinates. 
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Figure  125.  Profiles  of  Reynolds  Shear  Stresses  with  respect  to  Y <p 
X/L=0.8,  Re=4.0xl06,  a- 10°,  Local  Free-Stream  Coordinates 
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Comparison  of  u~  Reynolds  Stress  in  Regions  of  Converging 
Streamlines,  Re  =  4.0xl06,  a  =  10°,  DFVLR  Data  at  Re  =  7.2xl0^ 
(Kreplin  and  Vollmers,  1981). 
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Figure  127.  Comparison  of  v-  Reynolds  Stress  in  Regions  of  Converging 
Streamlines,  Re  =  4.0xl06,  DFVLR  Data  at  Re  =  7.2xl06  (Kreplin  and 
Vollmers,  1981). 
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Figure  128.  Comparison  of  w2  Reynolds  Stress  in  Regions  of  Converging 
Streamlines,  Re  =  4.0xl06,  a=I0°,  DFVLR  Data  at  Re  =  7.2xl06 
(Kreplin  and  Vollmers,  1981). 
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Figure  129.  Comparison  of  av  Reynolds  Stress  in  Regions  of  Converging 
Streamlines,  Re  =  4.0xl06,  a  =  10°,  DFVLR  Data  at  Re  =  7.2xlOft 
(Krepiin  and  Vollmers,  1981). 
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Figure  131 


Comparison  of  nv  Reynolds  Stress  in  Regions  of  Convergim 
Streamlines,  Re  =  4.0xl06,  a  =  10°,  DFVLR  Data  at  Re  =  7.2xl0 
(Kreplin  and  Vollmers.  1981). 
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Figure  133.  Profiles  of  Total  Turbulent  Shear  Stress  r,  X/L=0.8,  Re=4.0xl06, 
a=10°,  Local  Free-Stream  Coordinates. 
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Figure  134.  Profiles  of  Townsend’s  Structural  Parameter  Al.  X/L  =  0.7,  Re=4.0x  106, 
a-- 10°,  Local  Free-Stream  Coordinates. 
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Figure  135.  Profiles  of  Townsend’s  Structural  Parameter  A1 ,  X/L  =  0.8,  Re=4.0x  106, 
a=10°,  Local  Free-Stream  Coordinates. 
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Figure  137.  Profiles  of  X  Eddy  Viscosity  vx,  X/L=0.8,  Re=4.0x  106.  a=  10°,  Local 
Free-Stream  Coordinates. 
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Figure  139. 


Profiles  of  Z  Eddy  Viscosity  vz,  X/L=0.8,  Re=4.0x  106,  a  = 
Free-Stream  Coordinates. 


10°,  Local 


Illustrations 


1 Q  t 


/«/„*•) 


0.025 


0.020  r 


0.015  E 


0.010 


0.005 


top 


0.000 


■AB 


<p  =  20° 

9?  =  30° 
9?  =  40‘ 
(^  =  50c 


^  l  l 


0.1 


V6 


99 


Figure  140.  Profiles  of  vx/(Ure<5*),  X/L=0.7,  Re=4.0xl06,  a=10°,  Local  Free- 
Stream  Coordinates. 
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Figure  141.  Profiles  of  Vx/(Ure6*),  X/L=0.8,  Re  =  4.0xl(/\  a=10°,  Local  Free- 
Stream  Coordinates. 
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Figure  142.  Profiles  of  Anisotropy  Constant  N,  X/L  =  0.7,  Re  =  4.0xl06,  a=10°. 
Local  Free-Stream  Coordinates. 
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Figure  143.  Profiles  of  Anisotropy  Constant  N,  X/L  =  0.8,  Re  =  4.0xl06,  a=10°. 
Local  Free-Stream  Coordinates. 
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Figure  144.  Profiles  of  Mixing  Lengih  Lm/Lm:D,  X/L  =  0.7,  Re  =  4.0xl0<s.  a  =  10 
Local  Free-Stream  Coordinates. 
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Figure  145.  Profiles  of  Mixing  Length  Lm/Lm;o,  X/L=0.8,  Re  =  4.0xl06,  a=10°, 
Local  Free-Stream  Coordinates. 
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Figure  146.  Profiles  of  Flow  Gradient  Angle  ag.  X/L=0.7,  Re  =  4.0xl06.  q  =  10 
Local  Free-Stream  Coordinates. 
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Figure  147.  Profiles  of  Flow  Gradient  Angle  a£,  X/L  =  0.8,  Re  =  4.0xl06,  cr  =  10 
Local  Free-Stream  Coordinates. 
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Figure  148.  Profiles  of  Shear  Stress  Angle  as,  X/L  =  0.7,  Re  =  4.0x  106,  a=  10°,  Local 
Free-Stream  Coordinates. 
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Figure  149.  Profiles  of  Shear  Stress  Angle  as.  X/L  =  0.8.  Re  =  4,0x  106,  a=  10°,  Local 
Free-Stream  Coordinates. 
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Figure  150.  Profiles  of  Difference  Between  Flow  Gradient  Angie  and  Shear  Stress 
Angle  ag-as,  X/L=0.7,  Re=4.0xl06,  a=l0°,  Local  Free-Stream 
Coordinates. 
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Figure  151.  Profiles  of  Difference  Between  Flow  Gradient  Angle  and  Shear  Stress 
Angle  ag-a6,  X/L=0.8,  Re=4.0xl06,  a=10°,  Local  Free-Stream 
Coordinates. 
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Table  1. 


DFVLR,  West  Germany,  Measurements  and  Flow  Conditions 
6: 1  Prolate  Spheroid 


Measurements 

Reynolds  Number 

Angle  of  Attack  (deg) 

Skin  friction 

1.6x  106 

0, 

2.5,  5,  10.  30 

7.2X106 

Surface  oil-flow 

1.6xl06 

0. 

2.5.  5,  10,  30 

visualizations 

7.2  xlO6 

Surface  pressures 

1.6x  106 

0, 

2.5,  5,  10,  30 

7.2xl06 

Boundary  layer  profiles 

1.6x  106 

10.  30 

7.2X106 

Mean  velocity  field 

1.6X106 

10,  30 

7.2  xlO6 

Reynolds  stresses  in 

7.2xl06 

10 

boundary  layer 
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Table  2.  Boundary  Layer  Parameters,  X/L=0.7 


<p 

<$99  (cm) 

S'  (cm) 

0  (cm) 

H 

20° 

0.870 

0.121 

0.085 

1.424 

o 

0 

1.041 

0.159 

0.119 

1.336 

o 

O 

1.121 

0.174 

0.122 

1.426 

50° 

1.677 

0.315 

0.212 

1.486 

60° 

1.773 

0.304 

0.219 

1.388 

70° 

1.923 

0.246 

0.187 

1.316 

0 

O 

oc 

1.404 

0.164 

0.125 

1.312 

o 

o 

o 

1.263 

0.165 

0.120 

1.375 
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Table  3.  Boundary  Layer  Parameters,  X/L  =  0.8 


<p 

(cm) 

S'  (cm) 

0  (cm) 

H 

20° 

1.072 

0.194 

0.132 

1.470 

o 

0 

1.718 

0.341 

0.228 

1.496 

0 

O 

1.802 

0.422 

0.266 

1.586 

L/i 

o 

0 

2.258 

0.474 

0.315 

1.505 

60° 

2.470 

0.433 

0.311 

1.392 

o 

O 

r~ 

2.669 

0.384 

0.292 

1.315 

0 

O 

oo 

2.038 

0.235 

0.175 

1.343 

90° 

1.699 

0.234 

0.169 

1.385 
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Appendix  A.  Yawhead  Probe  Data  Reduction 


The  data  reduction  method  and  equations  for  the  mean  flow  measurements 
performed  using  the  three-dimensional  yawhead  probe  are  taken  from  Sung  (1985).  At 
each  data  point,  the  pressures  at  each  of  the  five  ports  on  the  probe,  shown  in  Figure  10. 
were  measured  and  recorded,  along  with  the  tunnel  free-stream  temperature  and  pressure. 
Once  the  data  were  taken  at  a  measurement  point,  the  following  sets  of  equations  were 
used  to  obtain  the  three  components  of  the  velocity  in  probe  coordinates. 

The  dimensionless  pressure  coefficients  are  given  by: 


Puxal 


{P*  -  P>)  __ 

A 

{Pi  ~  ^3) 

A 

(Pi  -  ^r)  _ 
A 

(Ps^  ~  Ps) 

A 


/,(“>  P) 
/2(«.  P) 
/3(a»P) 

=  /4(a.P) 


(A.  1) 
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where  P,  through  P5  refer  to  the  pressures  at  the  five  ports,  PT  is  the  total  pressure.  Ps 
is  the  static  pressure,  a  is  the  pitch  angle,  ft  is  the  yaw  angle,  and  A  and  Ps  are  given 
by: 


A  = 


P  -  P 

r\  rSw 


(P2  +  Pl  +  Pt  +  P5 


S^vg 


(A. 2) 


The  functions  f,  through  f4  are  calibration  curves  determined  by  Sung  (1985). 

After  the  pressures  P!  through  P5  are  measured,  the  quantities  A,  PSavg,  Cp  pilch. 
and  Cp  yav,  are  determined.  The  calibration  curves  f,  and  f:  are  used  to  determine  a  and 
ft.  Then  PT  and  Ps  are  found  using  the  calibration  curves  f3  and  f4.  The  magnitude  of 
the  velocity  vector  U[ot  is  found  using  the  incompressible  form  of  Bernoulli's  equation  as 
follows: 


U., 


(A. 3) 


The  components  of  the  velocity  vector  in  probe  coordinates  are  calculated  from: 

U  =  |  Uwt  |  cos  a  cos  P 

V  =  I U |  sin  P  (A. 4) 

W  =  I U a,  i  sin  a  cos  (1 
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Appendix  B.  X-Wire  Data  Reduction  Methods 


Appendix  B.l  X-Wire  Response  Equations 


This  appendix  describes  the  techniques  used  with  the  x-wire  hot-wire  anemometer 
to  obtain  mean  velocity  and  turbulence  quantities.  As  shown  in  Figure  6.  the  DANTEC 
55P61  probe  has  two  sensors,  denoted  by  subscripts  1  and  2,  arranged  nominally  at 
0=  ±45°  with  respect  to  the  axis  (probe  axis).  According  to  Jorgensen  (1971),  the 
relationship  between  the  effective  cooling  velocity  of  a  sensor  and  the  normal,  tangential, 
and  out-of-plane  cooling  velocities  is  given  by: 


ul*  *  kWB 


(B.l) 


The  tangential  and  normal  coefficients,  kT  and  kN.  were  initially  estimated  by  placing  the 
probe  at  different  orientations  in  a  TS1  Model  1 125  calibrator,  assuming  that  0=  ±45°. 
The  resulting  values  were  compared  with  those  obtained  by  Ha  (1989)  for  an  identical 
DANTEC  55P61  probe,  and  found  to  be  within  15 %  of  that  data.  Therefore,  because 
the  calibrations  from  Ha  were  made  using  a  more  precise  set-up  for  measuring  the  pitch 
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and  yaw  angles  of  the  probe,  the  values  from  Ha  (1989)  were  used  in  the  data  reduction. 
Because  the  coefficients  are  dependent  upon  the  orientation  of  the  sensor  with  respect  to 
the  free-stream  velocity  and  the  present  data  acquisition  and  reduction  method  utilizes 
rolling  of  the  probe  about  its  axis  to  obtain  all  six  Reynolds  stress  quantities,  curve  fits 
of  kT  and  kN  with  respect  to  yaw  and  pitch  angles  were  made  with  the  data  from  Ha  and 
used  during  the  data  reduction. 

Before  and  after  acquiring  data  for  each  run,  the  probe  was  calibrated  in  the 
Stability  Wind  Tunnel  by  placing  the  probe  axis  parallel  to  the  tunnel  free-stream 
velocity.  The  tunnel  velocity  was  run  at  different  settings  between  0  fps  and  180  fps,  and 
at  each  velocity  the  linearized  output  voltages  from  the  sensors  were  stored.  A  linear 
least  squares  curve  fit  was  made  for  each  calibration,  of  the  form: 

E  =  A  ~  BV  (B.2) 

where  E  is  the  output  voltage,  V  is  the  effective  cooling  velocity,  and  A  and  B  are 
constants.  In  each  case,  the  correlation  coefficient  was  greater  than  0.999.  The  tunnel 
temperature  and  static  pressure  were  also  stored  for  each  calibration. 

Tor  each  measurement  location,  the  probe  was  oriented  to  five  different  probe  roll 
angles  shown  in  Figure  6  and  mean  and  rms  voltages  for  each  sensor  were  recorded. 
After  acquiring  all  the  data  for  a  run,  the  voltages  were  converted  to  effective  cooling 
velocities  for  each  sensor  using  the  linearized  calibrations,  and  corrections  for  temperature 
changes  were  made  using  the  method  from  Larsen  and  Busch  (1974).  These  involved 
corrections  of  the  constant  B  in  equatio..  (B.2)  for  each  data  point,  of  the  form: 


Appendix  B 


211 


where  p  is  the  tunnel  static  pressure,  p0  is  the  calibration  reference  pressure,  n  is  the 
exponent  from  King’s  Law  determined  during  the  calibrations,  6,  is  a  parameter  equal 
to  1.09  according  to  Larsen  and  Busch,  and  e  is  given  by: 


e 


In  this  expression.  Tg  is  the  tunnel  free-stream  temperature,  T0  is  the  calibration  reference 
temperature,  and  Tw  is  the  temperature  of  the  hot-wire  sensors. 

The  resulting  velocities  consisted  of  the  mean  effective  cooling  velocity  U ^  and 
the  fluctuating  effective  cooling  velocity'  u\g  for  each  sensor  and  the  correlation  velocity 


ueffiuef2'  at  eac^1  Pr°be  roll  angle  ijr. 

In  order  to  obtain  all  the  mean  velocity  and  turbulence  measurements  at  a  point, 
equations  for  the  effective  cooling  values  must  be  solved.  The  following  derivation  gives 
the  required  equations. 

According  to  the  geometry  of  the  sensor  in  Figure  6  (Acrivlellis,  1978): 

UN  =  (U  +  u) cos0  *  [(K+  v)cosiJj  +  {W *  w)simjrlsin0 

UT  =  (U+u)sinQ  +  [(F+v)cosi{j  +  (IT+  w)sim|j]cos0  (B.3) 

UB  =  -(V  +  v)sim|r  ~  (W+  w)cost|f 


Now  define: 
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A  =  cos2©  +  kj  sin2© 

B  =  sin20  +  ^sm20 

n 

D  =  2cos0sin0  -  2£^cos0sin0 
=  (l  -  k})sin2Q 
F  =  sin2©  +  k2T cos2©  -  k“N 
-  B  -  k2N 


Substitute  equations  (B.3)  and  (B.4)  into  equation  (B.l)  to  obtain: 

Ua  -  (U*u)~ A 

+  (V ~  v)2|jBcos2t|f  +  k“N  sin2\{;| 

+  (W+  u’)2|5sin24r  +  k2N  sin2  41] 

+  (f/+u)(F+v)Dcosi|» 

-r  (U  +  u){V+\>)D sin4» 

+  (K-\’)(W/+H’)Fsin2t|r 


Expand  the  products  in  equation  (B.5),  and  assume  u2,  v2,  w2.  uv,  uw,  and  vw 
smaller  than  U.  Now  define  the  following: 

A'  =  A 

B  =  ficos'ty  +  kffSW \|» 

C'  =  5sin24»  +  ^,cos2vJf 
D'  =  Dcosij/ 

E'  =  DsiniJ/ 

F'  =  Fsin2ty 

and  substitute  equation  (B.6)  into  equation  (B.5)  to  obtain: 


(B.4) 


(B.5) 


are  much 


(B.6) 
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u]ff  =  A'V2  4  2A'Uu  +  B'F2  *  25'Fv 
+  C'W1  +  2C/PFw  +  D'C/F  4  D'C/v 

4  D  uF  4  £'£/jy  4  E'lJw  4  £'W1F 

4  F'vw  4  F/Fw  4  F'vW' 


(B.7) 


Factor  out  A'U2  in  equation  (B.7)  and  use  binomial  expansion  to  find  the  square  root, 
neelectine  hieher  order  terms.  The  result  is: 


U  =  jA'(U~u)  4  — (F* v)  4  — — -  (fF+  H’) 

2V^47  2vU; 

r  V  r  IF 

- —  ( F  4  2  \>)  4  — — — -  ( PF  +  2  H-) 


2v/4/  U 


2yAl  U 


— —  (FTF-  Fw  4  vBO - - — (D'V+E'W) 


2\j  A  '  U 


2\[aJ  U 


(B.8) 


Take  the  time  mean  and  neglect  smaller  terms  in  equation  (B.8)  to  obtain: 


U  „  =  <A'  U  -  V - —  W 


2JX1  2\A‘ 


(B.9) 


To  find  the  turbulence  quantity,  use  the  fact  that: 


Ueff  ~  Ueff 


(B. 10) 


to  obtain: 


«W  =  V"  + 


7  D'V-E'W 


2\[a/  U 


B V  F'W 


2  \jA '  sjA'U  2  vA7  U 


E1  C'W  F'V 

-  +  -  +  - 

2 {a*  {a'  U  2  JZ'U 


(B.ll) 
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Define  the  terms  in  brackets  in  equation  (B.  11)  as  P,  Q,  and  R  to  obtain  the  following 
for  each  sensor: 


Ueffl  =  P\u  +  <?!V  + 

“eff:  =  P2“  +  Q2v  + 


(B. 12) 


Square  each  of  the  expressions  in  equation  (B.12)  and  take  the  time  mean,  and  also  find 
the  product  of  the  two  expressions  in  equation  (B.12)  and  take  the  time  mean  to  obtain: 

=  P[u2  -  q\v2  +  R[w2  ^2PlQluv  +  2P:Rluw~2QlRlvw 
ueff2  =  P2u2  ~Q2v2  +  fcW2  +  2P2Q2uv  +  2P2R2uw  +  2Q1R2vw 

UefflUeff2  =  PXP2U  +Q\Q2V  +{P\Q2  +  P2@l)UV 

~^lR  +  P2Rl)uH^lplR2~Q2Rl)w 

Equation  (B.9)  contains  the  three  desired  mean  velocities  U,  V,  and  W.  To  solve 
for  these  quantities,  three  independent  equations  are  necessary.  The  three  equations 
(B.13)  contain  the  six  desired  Reynolds  stress  quantities.  In  order  to  solve  for  these 
quantities,  six  independent  equations  are  necessary'.  Therefore,  the  x-wire  probe  was 
rolled  to  the  five  different  angles  shown  in  Figure  6.  At  each  roll  angle,  as  previously 

described,  U ^  and  u 2rS  for  each  sensor  and  the  correlation  were  obtained. 

To  obtain  estimates  of  the  three  components  of  the  mean  velocities,  the 
expressions  for  U ^  (equation  B.9)  for  each  sensor  at  the  various  roll  angles  were 

grouped  as  follows: 

Group  1:  77^=0°),  77^(*  =  90°),  77^(i|r=  180°) 

Group  2:  77^(l(f=0o),  £/^(*=90°),  77^(i|r=  180°) 
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Group  3: 
Group  4: 


^(*=45°),  t/^(^  =  90°)-  ^(*  =  135°) 

I77^(^=450).  77^:(4r=90o),  T7^(4r=  135°) 

The  three  equations  in  each  of  the  four  groups  were  solved  simultaneously  to  find 
estimates  of  U,  V,  and  W,  and  these  estimates  were  averaged  to  obtain  the  final  values. 
The  estimates  for  each  mean  velocity  agreed  to  within  approximately  3%  of  the  averages. 

To  find  the  turbulence  quantities,  the  expressions  in  equation  CB.  13)  for  each 
sensor  at  the  various  roll  angles  were  grouped  as  follows: 


Group 

1: 

11 '81 

(4 '  = 

=  0°). 

ulnW=  45 

°). 

Urffl 

(r- 

=  90°), 

11 '81 

0|r  = 

=  135° 

)•  = 

18C 

>°), : 

Uf81 

=0°) 

Group 

T  • 

■ 

Uc82 

(♦  = 

=0°), 

°). 

u\m 

(♦s 

=90°), 

U'ff2 

(¥  = 

=  135° 

')«  U'82^y = 

180°). 

u '81 

=0°) 

Group 

3: 

U'8i 

W  = 

=0°), 

w^(4r=45 

°), 

u*ffl 

=  90°), 

11 '81 

(4>  = 

=  135° 

)•  = 

18C 

n. ' 

=45°) 

Group 

4: 

U'82 

(4r  = 

=0°), 

w^2(4'=45 

°). 

(4r  = 

-90°), 

U'ffZ 

(Vs 

=  135° 

18C 

»°), ; 

=45°) 

Group 

5: 

ulffl 

(4r  = 

=  0°), 

U'81  =43 

°). 

4 

(*  = 

=  90°), 

u m 

(¥  = 

=  135° 

).  Ueffl(\ Jr  = 

18C 

i°). : 

U'81 

=  90°) 

Group 

6: 

Ulff2 

(*  = 

=  0°), 

w^(^=45 

°). 

2 

U'8> 

(*  = 

=  90°), 

u'82 

(♦  = 

=  135° 

)>  ^(4^  = 

i8o°), ; 

u '81 

=90°) 

Group 

7: 

u'8l 

(t  = 

=0°), 

“tfj(¥=45 

°), 

“<81 

(*  = 

=  90°), 

u'81 

(♦  = 

=  135° 

)’  *4i(  +  = 

i8o°). : 

U'8l 

= 

=  135°) 

Group 

8: 

U'82 

(♦  = 

=0°), 

w^(4;=45 

°), 

(♦  = 

=  90°), 

u'82 

(♦  = 

=  135° 

)’  M^2(t  = 

i8o°), : 

U'82^Z 

=  135°) 

Group 

9: 

Ulffl 

(+  = 

=0°), 

“^(4»=45 

°), 

u*l 

(,|r  = 

=90°), 
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Group  10:  z^(iJj  =  0o),  r  =  45°),  =  90°), 

«^0|f  =  135°).  u2eff2{Tif=  180°),  =  180°) 

The  six  equations  in  each  group  were  solved  simultaneously  to  obtain  a  total  of  ten 
estimates  for  each  of  the  desired  turbulence  values,  and  these  estimates  were  averaged  to 

obtain  the  final  values.  The  estimates  for  u2 ,  v2,  and  w2  agreed  to  within  6%  of  the 
averages.  The  estimates  for  uv  and  uw  agreed  to  within  129c  of  the  averages,  and  the 

estimates  for  wv  agreed  to  within  20%  of  the  averages. 

After  the  desired  values  were  obtained  in  probe  coordinates,  they  were 

transformed  to  tunnel  coordinates  using  a  tensor  transformation,  taking  into  account  the 

pitch  angle  aH  of  the  probe  with  respect  to  the  XT  axis,  as  shown  in  Figure  6.  The 

equations  for  this  transformation  are  given  below,  in  which  the  subscript  H  refers  to  the 

hot-wire  probe  coordinates  and  the  subscript  T  refers  to  tunnel  coordinates: 

UT  =  UHcosaH  +  Fwsinafl 
VT  =  -U  H  sinaw  +  F^cosa^ 

*7  =  *H_  _  _ 

u2  t  =  u2  H  cos2aw  +  2uvwcosa  wsina  w  +  v2Hsin2aH 
v2 T  =  -  2uvh cosawsinaK  +  v2Hcos2aH 

=  w\,  _ 

u\’T  =  (v2w  -  u2 ^ cosa wsinaw  +  uvH(cos2aH  -  sin2aw) 
uwT  =  uwH  cosaff  +  vw^sina^ 
iWj.  =  -uw  H  sinaw  +  vw^cos  a  H 
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Appendix  B.2  Data  Reduction  FORTRAN  Program 


C  XW1RE.FOR  JANUARY  16,  1990. 

C 

C  THIS  PROGRAM  SOLVES  FOR  THE  SIX  REYNOLDS  STRESSES  AND  THE  THREE 
C  MEAN  VELOCITY  COMPONENTS  AT  A  POINT.  WITH  MEAN  AND  RMS  VOLTAGES 
C  AS  INPUTS.  IT  SOLVES  THE  EQUATIONS  SIMULTANEOUSLY  USING 
C  SUBROUTINE 

C  THE  INPUT  VOLTAGES  ARE  CORRECTED  FOR  TEMPERATURE  DIFFERENCES, 

C  USING  THE  METHOD  FROM  LARSEN  AND  BUSCH  (1974).  THE  QUANTITIES 
C  ARE  WRITTEN  TO  A  DATA  FILE  IN  PROBE  COORDINATES. 

C 

REAL  UEFF1(10).UEFF2(  10),UBAR(  10i.VBAR(10).WBAR(  10) 

REAL  UPEFF1(  10),UPEFF2(10),U  12PEFF(10),UP2(  10),VP2(  10),WP2(10) 

REAL  UVP(  10),UWP(10).VWP(10),PHI(5) 

REAL  A(2,2),B(2).X(2),AA(6,6),BB(6),XX(6) 

REAL  AP1(5).BP1(5).CP1(5),DP1(5).EP1(5),FP1(5) 

REAL  AP2(5).BP2(5).CP2(5),DP2(5).EP2(5),FP2(5) 

REAL  Pl(5).Ql(5).Ri(5).P2(5),Q2(5),R2(5) 

REAL  KN1(5).KN2(5).KT1(5).KT2(5).NCAL1,NCAL2 
REAL  A1(5),A2(5).B1(5),B2(5),D1(5),D2(5),F1(5),F2(5) 

INTEGER  ORDER(2),ORDERl(6) 

C 

C  BETA1  IS  GIVEN  IN  LARSEN  AND  BUSCH  (1974) 

C  TWIRE  IS  HOT-WIRE  TEMPERATURE 
C  PROMPTS  FOR  CONDITIONS  AT  CALIBRATION 
C 

BETA1  =  1.09 
TWIRE  =  200.0 

W'RITE(*,*)  ’ENTER  CALIBRATION  TEMP(F),  PRESSURE(PSI),  N1  AND  N2’ 
READ/*,*)  TCAL.PCAL,NCAL1,NCAL2 
TCAL  =  (TCAL-32.0)*5. 0/9.0 

READS  CALIBRATION  CONSTANTS  AND  LOCATIONS  OF  MEASUREMENT. 

READ(7,*,END  =  999)  AW1  ,BW1  .AW2.BW2 
WRITE!*,*)  ’ENTER  Y-LOCATION  OF  LOWEST  PT.,  IN  CM’ 

READ/*,*)  YFINAL 

WRITE(*,*)  ’ENTER  DISTANCE  FROM  LOWEST  PT.  TO  SURFACE,  IN  CM’ 
READ/*.*)  DIST 
DIST  =  YFINAL-DIST 

W'RITE/*,*)  ’ENTER  Z-LOCATION  IN  INCHES’ 

READ(*,*)  ZLOC 
PI  =  3. 141592654 
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INITIALIZE  TANGENTIAL  AND  NORMAL  COOLING  COEFFICIENTS. 

DO  2  1  =  1.5 
KT 1  <1 ) = 0. 3 
KT2(I)=0.3 
KN  1(I)=  1 . 14 
KN2(I)=  1. 14 

2  CONTINUE 

READ  WIRE  ANGLES 

WRITE!*.*)  ’INPUT  THETA1  AND  THETA2  IN  DEGREES' 

READ)*.*)  THETA1.THETA2 
THETA1  =  THETA  1*PI/ 180.0 
THET  A2  =  THET  A2*PI/ 1 80.0 
CT1  =COS(THETAl) 

ST1  =SIN(THETA1) 

CT2  =  COS(THETA2) 

ST2  =  SIN(THETA2) 

CALCULATE  COEFFICIENTS  FOR  DATA  REDUCTION  EQUATIONS 
DO  3  1=1.5 

A1(I)  =  CT1  *CT1  +KTUI)*KTHI)*ST1*ST1 
B1(I)  =  ST1*ST1  +  KT1(I)*KTHI)*CTT*CTI 
D1(I)  =  2.0*(1.0-KT1(I)*KT1(I))*CT1*ST1 
F1(I)  =  BUI)-KNUI)*KN1(1) 

A2(I )  =  CT2*CT2  +  KT2(I)*KT2(I)*ST2*ST2 
B2(I )  =  ST2*ST2  +  KT2(I)*KT2(I  )*CT2*CT2 
D2d)  =  2.0*(I.0-KT2(I)*KT2(I))*CT2*ST2 
F2fI)  =  B2(l)-KN2(I)*KN2(I) 

3  CONTINUE 

READ  IN  VOLTAGES  AND  CONVERT  TO  EFFECTIVE  VELOCITIES 
5  DO  10  1=1.5 

READ(7,*,END  =  999)  YLOC.PHI(I) 

READ(7,*,END  =  999)  E1.E2.E1P.E2P 
READ(7,*.END=999)  E1E2P.E1SQ.E2SQ.E1E2 
UEFF1(I)  =  (E1-AW1)/BW1 
UEFF2fI)  =  (E2-AW2)/BW2 
UPEFF1(I)  =  E1P/BW1/BW1 
UPEFF2(I)  =  E2P/BW2/BW2 
U12PEFF(I)  =  E1E2P/BW1/BW2 
PHIfI)  =  PHI(I)*PI/ 180.0 
10  CONTINUE 

R£AD(7,*.END  =  999)  YLOC.UINF.TINF.PINF 

C  CORRECTION  FOR  TEMPERATURE  USING  LARSEN  AND  BUSCH  EQUATION 
C  INCLUDES  CORRECTION  FACTORS  FOR  B1  AND  B2 
C 
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TlNF  =  (TlNF-32.0)*5. 0/9.0 
EPSILON  =  (TINF-TCAL)/(TWIRE-TCAU 
BCORR  l=(l  .0-BETA 1  “EPSILON  )**NCAL  1  “PINF/PCAL 
BCORR2=(l. 0-BETA1  “EPSILON  )**NCAL2*PINF/PCAL 
DO  15  1=  1,5 

UEFFl(I)  =  UEFFl(I)/BCORRl 
UEFF2(I)  =  UEFF2(I)/BCORR2 
U  PEFF 1  (I )  =  U  PEFF 1  (I  )/BCORR  1  /BCORR  1 
UPEFF2(I)  =  UPEFF2(I)/BCORRl/BCORR2 
U 1 2PEFFII )  =  U 1 2PEFF(I  l/BCORR  1  /BCORR2 
15  CONTINUE 

CALCULATE  MORE  COEFFICIENTS  FOR  DATA  REDUCTION  EQUATIONS 

DO  20  1  =  1.5 
CP  =  COS(  PHI(Il) 

SP  =  SIN(PHI(D) 

AP1(I)  =  A1(I) 

BPl(I)  =  Bim*CP“CP  +  KNl(I)*KNl(I)“SP*SP 
CP1(I)  =  BHI)*SP*SP  +  KN1(I)*KN1(I)*CP*CP 
DP1(I)  =  D1(I)*CP 
EPHI)  =  D1(I)*SP 
FPKI)  =  FHI)*SIN(2.0*PHI(I)) 

AP2(I)  =  A2(I) 

BP2(I )  =  B2H)*CP“CP  +  KN 1  (I  >“KN  1  fI)*SP*SP 
CP2(I)  =  B2(I)*SP*SP  +  KN  1(I)*KN  l(l)“CP*CP 
DP2(I)  =  D2(I)“CP 
EP2(I)  =  D2(I)*SP 
FP2n)  =  F2(I)»SIN(2.0*PHI(l)) 

20  CONTINUE 
C 

C  LOOP  TO  CALCULATE  THE  MEAN  VELOCITIES  USING  LUDCMP  AND  SOLVLU  TO 
C  SOLVE  SIMULTANEOUS  EQUATIONS. 

C 

DO  58  LOOP=  1,4 
N  =2 

NDIM  =  2 

A(1.D  =  SQRT(AP1(1)) 

A(1.2)  =  DPl(l)/2.0/SQRT(APl(l)) 

A(2.1)  =  SQRT(AP2(1)) 

A(2.2)  =  DP2(  1  )/2.0/SQRT(  AP2(  1 )) 

B(1)  =  UEFF1(1) 

B(2)  =  UEFF2(1) 

CALL  LUDCMP( A, N, NDIM, ORDER) 

CALL  SOLVLU(  A, B,X,N, NDIM, ORDER) 

UBARO)  =  X(l) 

VBAR(1)  =  X(2) 

A(1.1)  =  SQRT(AP1(3)) 

A(1.2)  =  EPl(3)/2.0/SQRT(APl(3)) 

A(2, 1 ) =SQRT(  AP2(3)) 

A(2,2)  =  EP2(3)/2.0/SQRT(AP2(3)) 
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B(1)  =  UEFF1(3) 

B(2)  =  UEFF2(3) 

CALL  LUDCMP(A,N.NDIM, ORDER) 
CALL  SOLVLU(A.B.X.N.NDIM.ORDER) 
UBAR(2)  =  X(1) 

WBAR(1)  =  X(2) 

A(L1)  =  SQRT(AP1(5)) 

A(  1 .2)  =  DP  1  (5)/2.0/SQRT(  AP1(5)) 

A(2,l  )  =  3QRT(AP2f5)) 

At2,2)  =  DP2(5)/2.0/SQRT(  AP2(5)) 

B(1)  =  UEFF1(5) 

B(2)  =  UEFF2(5) 

CALL  LUDCMP(A.N.NDIM. ORDER) 
CALL  SOL VLU(A.B,X.N.NDIM. ORDER) 
UBAR(3)  =  X(1) 

VBAR(2)  =  X(2) 

A(  1 . 1)  =  SQRT(AP1(2)) 

A(  1 ,2)  =  DPl(2)/2.0/SQRT(APl(2)) 

A(2,l  )  =  SQRT(AP2(2)) 

A(2.2)  =  DP2(2)/2.0/SQRT(AP2(2)) 

B(1)  =  UEFF1(2) 

B(2)  =  UEFF2(2) 

C ALL  L L DCM P( A . N , NDLM .ORDER ) 
CALL  SOLVLU(A,B,X.N,NDIM.ORDER) 
UBAR(4)  =  X(1) 

VPW  =  X(2) 

A(1,1)  =  SQRT(APH4)) 

Af  1 ,2)  =  DPU4)/2.0/SQRT(APl(4)) 

AC2. 1 )  =  SQRT(  AP2(4)) 

A('2,2)  =  DP2(4)/2.0/SQRT(AP2(4)) 

B(1)  =  UEFF1(4) 

B(2)  =  UEFF2(4) 

CALL  LUDCMP(A.N.NDIM.ORDER) 
CALL  SOLVLU(A,B.X.N.NDIM, ORDER) 
UBAR(5)  =  X(  1) 

VMW  =  Xf2) 

VBAR(3)  =  (VPW  +  VMW)/2.0 

WB  AR(  2 )  =  ( VPW- VM  W)/2.0 

USUM  =0.0 

VSUM  =0.0 

WSUM  =0.0 

DO  30  1=1,5 

WR1TE(*.*)  ’Uest  =  ’.UBAR(I) 

USUM  =  USUM  +  UBAR(I) 

30  CONTINUE 
UAVG  =  USUM/5.0 
WR1TE(*,*)  ’Uavg  =  ’,UAVG 
DO  40  1=1,3 

WRITE(*,*)  ’Vest  =  \VBAR(I) 

VSUM  =  VSUM  +  VBAR(I) 

40  CONTINUE 
VAVG  =  VSUM/3.0 
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WRITE(*,*)  'Vavg  =  WAVG 
DO  50  [  =  1.2 

WRITE(*,*)  ’West  =  \WBAR(I) 

WSUM  =  WSUM  +  WBAR(l) 

50  CONTINUE 
W'AVG  =  WSUM/2.0 
WRITE!*.*)  ’Wavg  =  ’.WAVG 
C 

C  FIND  NEW  COOLING  COEFFICIENTS  BASED  ON  MEAN  VELOCITIES,  USING  CURVE 
C  FITS  TO  DATA  FROM  HA  (1989). 

C 

DO  55  1=1,5 
CP  =  COS(PHl(I)) 

SP  =  SIN(PHI(I)) 

UN  =  CT1*UAVG  +  ST1*CP*VAVG  +  ST1*SP*WAVG 
UT  =  -STl*UAVG-‘-CTl*CP*VAVG  +  CTl*SP*WAVG 

ub  =  -sp*vavg+cp*wavg 

YAW  =  ABS(ATAN(UT/UN)) 

PITCH  =  ABS(  AT  AN(UB/UN» 

K.Tl(I)  = -0.2572 1*YAW**3  +  1. 109*  YAW*  YAW- 1.691 16*YAW  +  1 . 15076 
KN  HI)  =-0.925408*PITCH**5  +  4.97674*PITCH**4-10.2698*PITCH**3  +  10.2 
1 2954  *  PITCH  *PITCH-5. 13504*  PITCH +  2. 19519 
UN  =  CT2*UAVG  +  ST2*CP*VAVG  +  ST2*SP*WAVG 
UT  =  -ST2*UAVG  +  CT2*CP*VAVG  +  CT2*SP*WAVG 
UB  =  -SP*  V  AVG  +  CP*  WAVG 
YAW  =  ABS(ATAN(UT/UN)) 

PITCH  =  ABS(  AT  AN(UB/UN)) 

KT2(I)  =  -0.2572 1  *YAW**3  +  1 . 109*  YAW*YAW-1 .691 16*YAW  +  1 .15076 
KN2(I)  =  -0.925408*PITCH**5  +  4.97674*PITCH**4-10.2698*PITCH**3  +  10.2 
1 2954*  PITCH ‘PITCH-5. 13504*PITCH +  2. 19519 


RECALCULATE  DATA  REDUCTION  COEFFICIENTS  • 

A1(I)  =  CT1*CT1  +  KT1(I)*KT1(I)*ST1*STI 
B 1  (I>  =  ST  1  *ST  1  +KT1(I)*KT1(I)*CT1*CT1 
D1(I)  =  2.0*(1.0-KTI(I)*KT1(I))*CT1*ST1 
F 1  (I)  =  B 1  (I  )-KN  1  (I)*KN  1  (I) 

A2(I)  =  CT2*CT2  +  KT2(I)*KT2(I)*ST2*ST2  • 

B2(I)  =  ST2*ST2  +  KT2(I)*KT2(I)*CT2*CT2 
D2(I)  =  2.0*(  1 .0-KT2(I)*KT2(I))*CT2*ST2 
F2(I)  =  B2(I  i-KN2(I)*KN2(I) 

AP1(I)  =  AI(I) 

BP1(I)  =  BI(I)*CP*CP+KN1(I)*KN  1(I)*SP*SP 

CP1(I)  =  B1(I)*SP*SP  +  KN1(I)*KN1(I)*CP*CP  • 

DP1(I)  =  D1(I)*CP 

EP1(I)  =  D1(I)*SP 

FP1(I)  =  F1(I)*SIN(2.0*PHI(I)) 

AP2(I)  =  A2(I) 

BP2(I)  =  B2(I)*CP*CP  +  KN1(I)*KN1(I)*SP*SP  _ 

CP2(I)  =  B2(I)*SP*SP+KN1(I)*KN1(I)*CP*CP  • 

DP2(I)  =  D2(I)*CP 
EP2(l)  =  D2(I)*SP 
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FP2(I)  =  F2(I)*SIN(2.0*PHI(I)) 

55  CONTINUE 

WRITEC.*)  ’END  OF  LOOP', LOOP 
58  CONTINUE 

CALCULATE  P.Q.R  (COEFFICIENTS  FOR  TURBULENCE  EQUATIONS) 

DO  60  1=1,5 
SQRAP1  =SQRT(AP1(I)) 

PHI)  =  SQRAP1+(DP1(I)*VAVG  +  EPHI)*WAVG)/2.0'SQRAP1/UAVG 

Ql(I)  =  DPl(l)/2.0/SQRAPl+BPl(I)*VAVG/SQRAPl/UAVG  +  FPKI)*WAVG/2.0/SQ 

1RAP1/UAVG 

Rl(I)  =  EPl(I)/2.0/SQRAPl  +  CPl(I)*WAVG/SQRAPl/UAVG  +  FPl(I)*VAVG/2.0/SQ 

1RAP1/UAVG 

SQRAP2  =  SQRT(AP2(I)) 

P2m  =  SQRAP2  +  (DP2(I)*VAVG  +  EP2(I)*WAVG)/2.0/SQRAP2/UAVG 

Q2(I )  =  DP2(I  J/2.0/SQRAP2  +  BP2(  I)*  V  A  VG/SQRAP2/U  A  VG  +  FP2(I  )*W  AVG/2.0/SQ 

1RAP2/UAVG 

R2(I)  =  EP2n)/2.0/SQRAP2  +  CP2(I)*WAVG/SQRAP2/UAVG  +  FP2(I)*VAVG/2.0/SQ 
1RAP2/UAVG 
60  CONTINUE 
ICOUNT  =  5 

THIS  LOOP  FINDS  THE  TURBULENCE  VALUES  FROM  VALUES  FOR  WIRE  SENSOR  2. 
A  SIMILAR  LOOP  CAN  BE  ADDED  FOR  SENSOR  1. 

DO  65  1=1,5 

AA(  1,1)  =  P2(  1  )*P2(  l ) 

AA(  1,2)  =  Q2(  1  )*Q2(  1 ) 

AA(  1,3)  =  R2(  1  )*R2(  1 ) 

AA(  1 ,4)  =  2.0*P2(  1  )*Q2(  1 ) 

AA(1,5)  =  2.0*P2(1)*R2(1) 

AA(1,6)  =  2.0*Q2(1)*R2(1) 

AA(2, 1 )  =  P2(2)*P2(2) 

AAf2,2)  =  Q2(2)*Q2(2) 

AA(2,3)  =  R2(2)*R2(2) 

AA(2,4)  =  2.0*P2(2)*Q2(2) 

AA(2,5)  =  2.0*P2(2)*R2(2) 

AA(2,6)  =  2.0*Q2(2)*R2(2) 

AA(3,l)  =  P2(3)*P2(3) 

AA(3,2)  =  Q2(3)*Q2(3) 

AA(3.3)  =  R2(3)*R2(3) 

AAf3,4)  =  2.0*P2(3)*Q2(3) 

AA(3,5)  =  2.0*P2(3)*R2(3) 

AA(3,6)  =  2.0*Q2(3)*R2(3) 

AA(4, 1  )  =  P2(4)*P2(4) 

AA(4,2)  =  Q2(4)*Q2(4) 

AA(4,3)  =  R2(4)*R2(4) 

AA(4,4)  =  2.0*P2(4)*Q2(4) 

AA(4,5)  =  2.0*P2(4)*R2(4) 

AA(4,6)  =  2.0*Q2(4)*R2(4) 

AA(5,1)  =  P2(5)*P2(5) 
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AA(5,2)  =  Q2(5)*Q2(5) 

AA<5.3)  =  R2(5)*R2(5) 

AA!5.4)  =  2.0*P2!5)*Q2(5) 

AA(5.5)  =  2.0*P2!5)*R2!5) 

AA(5,6)  =  2.0*Q2(5)*R2(5) 

AA(6,1)  =  P1(I)*P2(I) 

AA(6,2)  =  Q1(I)*Q2(I) 

AA(6.3)  =  R1(I)*R2<I) 

AA(6,4)  =  P1(1)*Q2(I)  +  P2(I)*Q1(I) 

AA<6.5»  =  P1!I)*R2(I)  +  P2(I)*R1(1) 

AA(6,6)  =  Q1(I)*R2!I)  +  Q2(I)*R1(I) 

BB(1)  =  UPEFF2(1) 

BB(2)  =  UPEFF2(2) 

BB!3)  =  UPEFF2(3) 

BB(4)  =  UPEFF2(4) 

BB(5)  =  UPEFF2(5) 

BB(6)  =  U12PEFF(I) 

N  =  6 
NDIM  =6 

CALL  LU'DCMPi  AA.N.NDIM. ORDER  1) 

CALL  SOLVLU(AA.BB.XX.N.NDIM.ORDERl) 

IF(XX(  1  ).LE.O.O.OR.XX(2).LE.O.O.OR.XX(3).LE.O.O)THEN 
ICOUNT  =  ICOUNT-1 
XX(1)=0.0 
XX(2)  =  0.0 
XX!3)=0.0 
XX(4)  =  0.0 
XX!5)  =  0.0 
XX(6)=0.0 

WRITE!*,*)  ’  NEGATIVE  TURBULENCE  *  *  *  TOSSED  OUT!!!!’ 

ENDIF 

UP2(I)  =  XX(1) 

VP2(Il  =  XX(2) 

WP2(I)  =  XX(3) 

UVP(I)  =  XX(4) 

UWPd)  =  XX(5) 

VWP(I)  =  XX!6) 

WRITE!*,*)  ’upT  =  *,UP2(I) 

WRITE(*,*)  ’vpA2  =  \VP2(I) 

WRITE!*,*)  ’wp~2  =  ’,WP2(I) 

WRITE!*,*)  ’upvp*2  =  ’.UVP(I) 

WRITE!*.*)  ’upwp~2  =  ’,UWP!I) 

WRITE!*,*)  vpwp*2  =  ’.VWP(I) 

65  CONTINUE 
C 

C  THE  ESTIMATES  ARE  AVERAGED  AND  WRITTEN  TO  THE  OUTPUT  DATA  FILE 
C 

WRITE!*,*) 

WRITE!*.*) 

UP2SUM  =0.0 
VP2SUM  =0.0 
WP2SUM  =0.0 
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L'VPSUM  =0.0 
UWPSUM  =0.0 
VWPSUM  =0.0 
DO  70  1=1.5 

UP2SUM  =  UP2SUN1  +  UP2(I) 

VP2SUM  =  VP2SUM  +  VP2(I) 

WP2SU  M  =  WP2SU  M  +  WP2(  I ) 

L’VPSUM  =  UVPSUM  +  UVP(I) 

UWPSUM  =  UWPSUM  +  UWP(I) 

V  W  PS  U  M  =  VWPS  U  M  +  V  WP(  I  j 
70  CONTINUE 
X  COUNT  =  FLOAT)  1COUNT) 
up:  AVG  =  UP2SUM /XCOUNT 
VP2AVG  =  VP2SUM/XCOUNT 
wp:  AVG  =  WP2SU  M  /X  COUNT 
UVPAVG  =  UVPSUM/XCOUNT 
UWPAVG  =  UWPSUM/XCOUNT 
VWP  AVG  =  VWPSUM/XCOUNT 
WRITE!*,*)  ’upA2avg  =  \UP2AVG 
WRITE!*.*)  ’vpA2avg  =  \VP2AVG 
WRITE!*,*)  ’wp~2avg  =  ’.WP2AVG 
WRITE(*.*)  'upvpavg  =  '.UVPAVG 
WRITE!*.*)  ’upwpavg  =  ’.UWPAVG 
WRITE!*,*)  ’vpwpavg  =  ’.VWPaVG 
YLOC  =  (YLOC-DIST)/2.54 
ALPHA  =  ATAN(V  AVG/UAVG) 

BET  A  =  AT  AN(W  AVG/U  AVG) 

ALPHAD  =  ALPHA*  180.0/PI 
BET  AD  =  BETA*  1 80.0/PI 

UR  =  SQRT(UAVG*U  AVG  +  V  AVG*VAVG  +  WAVG*WAVG) 

WRITE! 8.*  I  YLOC.ZLOC.UINF 
WRITEI8.*)  UR.ALPHAD.BETAD 
WRITE/8.*)  UAVG.VAVG.WAVG 
WRITEI8.*)  UP2AVG, VP2AVG.WP2AVG 
WRITE!  8 .  * )  UVPAVG,  UWPAVG  .VWP  AVG 
GOTO  5 
999  STOP 
END 
C 

C  THE  FOLLOWING  SUBROUTINES  SOLVE  SIMULTANEOUS  EQUATIONS. 

C  THEY  ARE  TAKEN  FROM  "APPLIED  NUMERICAL  ANALYSIS"  BY  CURTIS  GERALD 
C  AND  PATRICK  WHEATLEY. 

C 

SUBROUTINE  LUDCMP(A,N,NDIM, ORDER) 

REAL  A(NDIM.N) 

INTEGER  N.NDIM.ORDER(N) 

.  REAL  SUM 

INTEGER  I.KCOL.NM  1  ,JM  1  .JCOL.JPl ,IROW 
DO  10  1=  !,N 
OPDER(I)  =  I 
10  CONTINUE 

CALL  APVT(A,N,ND1M. ORDER.  1) 
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IF(ABS(A(1,1)).LT. l.OE-5)  THEN 
PRINT  100 
RETURN 
ENDIF 

DO  20  KCOL  =  2,N 

A(  1  ,KCOL)  =  A(  1  ,KCOL)/A(  1,1) 

20  CONTINUE 
NM1  =N-1 

DO  80  JCOL  =  2,NMl 
JM 1  =  JCOL-1 
DO  50  1R0W  =  JC0L,N 
SUM  =0.0 

DO  40  KCOL=  1.JM1 

SUM  =  SUM  +  At  lROW.KCOL)*A(KCOL.JCOL) 

40  CONTINUE 

A(IROW.JCOL)  =  A<IROW,JCOL)-SUM 
50  CONTINUE 

CALL  APVT(A.N.NDIM.ORDER.JCOL) 

IF( ABS(  A(JCOL,JCOL)).LT.  1 .0E-5)  THEN 
PRINT  100 
RETURN 
ENDIF 

JP1  =  JCOL+  1 
DO  70  KCOL  =  JP1  ,N 
SUM  =0.0 

DO  60  IROW  =  1  ,JM  1 

SUM  =  SUM  +  A(JCOL.IROW)*  A(IROW.KCOL) 

60  CONTINUE 

A(JCOL.KCOL)  =  ( ADC  '  L,KCOL)-SUM)/A(JCOL.JCOL) 

70  CONTINUE 
80  CONTINUE 
SUM  =0.0 

DO  90  KCOL=  l.NMl 
SUM  =  SUM  +  A(N,KCOL)*A(KCOL,N) 

90  CONTINUE 

AfN,N)  =  A(N,N)-SUM 
RETURN 

100  FORMAT!/. '  VERY  SMALL  PIVOT  ELEMENT  INDICATES  A  NEARLY 
+  ’  SINGULAR  MATRIX.  ’) 

END 


SUBROUTINE  APVT(A,N,NDIM, ORDER, JCOL) 
REAL  A(NDIM.N) 

INTEGER  ORDER(N),N,NDIM,JCOL 
REAL  SAVE, ANEXT, BIG 
INTEGER  IPVT.JP1  ,IROW,KCOL,ISAVE 
IPVT  =  JCOL 

BIG  =  ABS(  AfJCOL.JCOL)) 

IP1  =  JCOL+  1 
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DO  10  1ROW  =  JP1  ,N 

ANEXT  =  ABS(A(IROW,JCOL)) 
IF(ANEXT.GT.BIG)  IPVT  =  IROW 
10  CONTINUE 
IF(IPVT.EQ.JCOL)  THEN 
RETURN 
ENDIF 

DO  20  KCOL=  l.N 
SAVE  =  A(JCOL.KCOL) 

A(JCOL.KCOL)  =  A(IPVT.KCOL) 

A(IPVT.KCOL)  =  SAVE 
20  CONTINUE 

1SAVE  =  0RDER(JC0L) 

ORDER(JCOL)  =  ORDER(IPVT) 

ORDER1IPVT )  =  ISAVE 
RETURN 
END 
C 
C 
C 
C 
C 

SUBROUTINE  SOLVLU(LU,B,X,N.NDIM, ORDER) 
REAL  LU(NDIM,N),B(N),X(N) 

INTEGER  ORDER(N),N,NDIM 
REAL  SUM 

INTEGER  UJCOLTROW.NPl .NVBL.IM  1 
DO  10  1=1. N 
J  =ORDER(I) 

X(I)  =  B(J) 

10  CONTINUE 
X(1)  =  X(1)/LU(1.1) 

DO  50  IROW  =  2.N 
IM 1  =IROW-l 
SUM  =0.0 

DO  40  JCOL=  1.IM1 

SUM  =  SUM  +  LU(IROW\JCOL)*X(JCOL) 

40  CONTINUE 

XHRO  W)  =  (X(IROW)-SUM)/LU(IROW,IROW) 

50  CONTINUE 
DO  70  IROW  =  2,N 
N  VBL  =  N  -IRO  W  +  1 
SUM=0.0 
NP1  =NVBL+ 1 
DO  60  JCOL  =  NPl,N 
SUM  =  SUM  +  LU(NVBL,JCOL)*X(JCOL) 

60  CONTINUE 

X(NVBL)  =  X(NVBL)-SUM 
70  CONTINUE 
RETURN 
END 
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